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Abstract 



We have obtained nonperturbative one-loop expressions for the mean energy-momentum tensor 
and current density of Dirac's field on a constant electric-like background. One of the goals of 
this calculation is to give a consistent description of back-reaction in such a theory. Two cases 
of initial states are considered: the vacuum state and the thermal equilibrium state. First, we 
perform calculations for the vacuum initial state. In the obtained expressions, we separate the 
contributions due to particle creation and vacuum polarization. The latter contributions are related 
to the Heisenberg-Euler Lagrangian. Then, we study the case of the thermal initial state. Here, 
we separate the contributions due to particle creation, vacuum polarization, and the contributions 
due to the work of the external field on the particles at the initial state. All these contributions are 
studied in detail, in different regimes of weak and strong fields and low and high temperatures. The 
obtained results allow us to establish restrictions on the electric field and its duration under which 
QED with a strong constant electric field is consistent. Under such restrictions, one can neglect the 
back-reaction of particles created by the electric field. Some of the obtained results generalize the 
calculations of Heisenberg-Euler for energy density to the case of arbitrary strong electric fields. 

PACS numbers: 12.20.Ds,11.15.Tk,11.10.Wx 

1 Introduction 

It is well-known that quantum field theory (QFT) in external backgrounds is a good model for the study 
of quantum processes in cases when a part of a quantized field is strong enough to be treated as a 
classical one. Here, it is assumed that quantum processes under consideration do not affect significantly 
the classical field (back-reaction is supposed to be small) and the field is treated as an external one. Of 
course, it is well-understood that, in principle, back-reaction must be calculated, at any rate, to evaluate 
the limits of applicability of the obtained results. However, the latter problem often remains open due 
to its high complexity. From physical considerations, it is clear that back-reaction may be very strong, 
so it must be taken into account for external backgrounds that can create particles from vacuum (for 
theories with an unstable vacuum under the action of an external field) and produce actual work on 
particles. The effect of particle-creation from vacuum by an external field (vacuum instability in external 
fields) ranks among the most intriguing nonlinear phenomena in quantum theory. Its consideration is 
theoretically important, since it demands that one should go beyond the scope of perturbation theory, and 
its experimental observation would be able to verify the validity of a theory in the domain of superstrong 
fields. Concerning the history of the subject, the reader may refer to the books [TJ [51 [31 H] and to the 
articles [5 for recent computational developments. Typical problems examined in the models with an 
unstable vacuum are related to the calculation of the density of particles created from vacuum. In some 
cases, this allows one to make phenomenological conclusions about back-reaction; see, e.g., [B]. However, 
a complete description of back-reaction is related to the calculation of mean values for the current density 
and energy-momentum tensor (EMT) of the matter field. This problem, is, in fact, one of the goals of 
the present article. In this connection, we recall that Heisenberg and Euler calculated the mean value of 
energy density of Dirac's field in a constant electromagnetic background at zero temperature under the 
condition of weakness of the electric fielcQ E, |E| « £ c = M 2 c 3 /eh ~ 1,3 • 10 16 V/cm, see [7]. Having 
this energy density at hand, they constructed a one-loop effective Lagrangian C. It turns out that this 
Lagrangian had been re- introduced by Schwinger [5] (for subsequent developments, see the review [5]) to 
define the vacuum-to-vacuum transition amplitude c v as follows: 
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without any restrictions on the external field. Schwinger demonstrated that the probability P v of a 
vacuum to remain a vacuum in a constant electric field is related to the imaginary part of C as follows: 

P v = \ Cv \ 2 = exp{-VT2Im£} . (2) 

Here, T is the field duration, and V is the volume of observation. Obviously, the effect can actually be 
observed as soon as the strength of the external field approaches the critical value E c . Given this, there 
remains an open question of what is the form of the energy density, as well as the form of the entire EMT, 
in arbitrary constant fields. In the present article, we study the mean current density and EMT of a 
spinor field subject to the action of an electric-like external background, while paying a special attention 
to the initial state of thermal equilibrium. A more exact and detailed setting of the problem is given in 
the next section. 

Aside from the principal importance of describing back-reaction for an understanding of the limits of 
applicability of QED in external fields, the calculations of the article are closely related to a wide class of 
problems (in the framework of quantum-field models in external fields) that recently attract considerable 
attention. Below, we present a brief review of these problems and of their importance, which may be 
regarded as an additional motivation for the present study. 

Even though an actual possibility of creating superstrong fields under laboratory conditions does 
not exist at present, the e + e~-pair production from vacuum by a slowly varying external electric field is 
relevant to phenomenology with the advent of a new laser technology capable of accessing the true strong- 
field domain. It is widely discussed [10] at the SLAC and TESLA X-ray laser facilities. Such strong fields 
may be relevant in astrophysics, where characteristic values of electromagnetic and gravitational fields 
near black holes are enormous. The Coulomb barrier at the quark surface of a hot strange star may be 
a powerful source of e + e~-pairs, which are created in extremely strong constant electric fields (dozens of 
times higher that the critical field E c ) of the barrier and the flow away from the star [TT]. Such emission 
could be the main observational signature of quark stars. The electric field near a cosmic string can also 
be extremely strong [12] . 

Apart from pure QED problems, there exist some closely related QFT problems in which vacuum 
instability in various electric-like external backgrounds plays an important role, for example, phase tran- 
sitions in non-Abelian theories, the problem of boundary conditions, or the influence of topology on 
vacuum, the problem of a consistent vacuum construction in QCD and GUT, multiple particle-creation 
in the context of heavy-ion collisions, etc. Recently, it has also been recognized that the presence of a 
background electric field must be taken into account in string theory constructions; see, e.g., [13j and 
references therein. 

In the early 1970s, finite-temperature effects were recognized as very important in QFT, and then 
Dittrich [14] computed the one-loop effective potential of finite-temperature QED in the presence of a 
constant magnetic field, starting from the Schwinger relation ([I}. This was followed by an intense study of 
finite-temperature and finitc-chcmical-potcntial one-loop effects for fermions in a constant uniform mag- 
netic field; see for a review |15[ 116] . This study was motivated by its possible applications to various areas 
of physics, including astrophysics, condensed matter, particle physics, etc., where strong electromagnetic 
fields may exist. The presence of an electric-like field implies a qualitative change (in comparison with a 
magnetic-like external field alone) in the character of processes involved, due to vacuum instability and 
the work of the field on the particles at the initial state. 

In some of the above-mentioned examples with an electric-like external field, the initial state is in 
thermal equilibrium. In the chromoelectric flux-tube model [17] , the back-reaction of created pairs in- 
duces a gluon mean field and plasma oscillations (see [6] and references therein). Various problems of 
cosmological QCD phase-transitions and dark matter formation arc discussed on the basis of the chromo- 
electric flux tube model (see, for example, [18j and references therein). It appears that for a calculation 
of particle creation in this model one needs to apply the general formalism of QFT for pair-production 
at finite temperatures and at zero temperature, both from vacuum and from many-particle states (see 
the corresponding physical reasons in [19l [20] [21] ) . A consideration of various time scales in heavy-ion 
collisions shows that the stabilization time (a time interval during which the characteristic asymptotic 
form for differential mean numbers of particles created by a constant chromoelectric field from vacuum 
is achieved) is far less than the period of plasma and mean-field oscillations. Then, the approximation 
of a strong quasiconstant chromoelectric field can be used in the treatment of such collisions during a 
period when the partons produced can be regarded as weakly coupled, due to the property of asymptotic 
freedom in QCD. It may also be reasonable to neglect dynamical back-reaction effects [22]. In the case 
of a strange star [11] , it has been argued that there is a macroscopic time interval during which local 
thermal equilibrium is achieved, but the strong electric field has not yet been depleted, due to e + e _ -pair 
production. By this reason, there has been a considerable interest in establishing a relevant formalism 
for finite-temperature QED with an unstable vacuum. An adequate techniques for such a case has been 
proposed in [23] . The special case of thermally-influenced pair production in a constant electric field has 
been studied at the one-loop level in [HI [22 [H 123 [26] . The results of [H [H [25] , obtained by various 
methods in the framework of the generalized Furry representation, are in mutual agreement. The authors 
of |I9[ 126] have arrived at different results, being, at the same time, in contradiction with themselves. 
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We believe that the difficulties [TH1 [53] are related to some attempts to generalize relation ([T]) to nonzero 
temperatures; see Discussion and Summary. 

The paper is organized as follows. In Section 2, we describe the setting of the problem and recall 
some necessary details of the non-perturbative techniques that we use in our calculations. In Section 3, 
we calculate the current density and EMT for the vacuum initial state. In the obtained expressions, we 
separate the contributions due to particle creation and vacuum polarization. The latter contributions 
are related to the Heisenberg-Euler Lagrangian. In Section 4, we calculate the current density and EMT 
for the thermal initial state. In the obtained expressions, we separate the contributions due to particle- 
creation, vacuum polarization, and the contributions that appear due to the work of the external field on 
the particles at the initial state. All these contributions are studied in detail in different regimes, limits 
of weak and strong fields and low and high temperatures. We have established restrictions on the electric 
field and its duration under which QED with a strong constant electric field is consistent. Within these 
restrictions, one can neglect the back-reaction of particles created by the electric field. In the last section, 
we summarize and discuss the main results. Some necessary but cumbersome calculations are placed in 
the Appendix. 

In what follows, we use the relativistic units h = c = 1 in which the fine structure constant is 
a = e 2 /ch = e 2 . 

2 The setting of the problem 

We consider a quantized spinor (Dirac) field ip(x), x — (a: = £, x), in an external electromagnetic back- 
ground specified by a quasiconstant electric field and a constant magnetic field. The field strength F^ v 
has nonzero invariants. In this case, we can always choose a reference frame in which the electric and 
magnetic fields are parallel and directed along the a; 3 -axis: 

F„ = F^ + F^, F^=E(Sl5l-8 3 X), F* v = B^l - Sffi) . (3) 

Our aim is to study the mean values of electric current and EMT in states that have evolved from 
vacuum and thermal equilibrium. It is well-known that a uniform constant electric field acting during 
an infinite time creates an infinite number density of pairs from vacuum and produces infinite work on 
particles. There arise other divergences in the course of QED calculations related to this phenomenon. 
This is why one needs a regularization to deal with these divergences. One type of such a regularization 
is the following, namely, we consider, from the beginning, a uniform electric field, which efficiently acts 
only during a sufficiently large but finite time T. Such a regularization has been used in our works 
[21 [UJ I n these works, we studied particle-creation effects in the so-called T-constant field, being a 
uniform electric field, constant (and nonzero) within a time interval T, and zero outside this intervaQ- 
A similar regularization is used in the present article: we choose the T-constant field to switch on at 
the instant t\ — —T/2 and to switch off at the instant t 2 = T/2. The results of our calculations will 
be presented in a form in which contributions of different nature are separated. These are contributions 
due to vacuum polarization, particle-creation from vacuum, and the contributions due to the work of the 
external field on the particles at the initial state. Some of these contributions have constant parts, as 
well as parts depending on the time interval t — t\ that passes since the instant the electric field turns on. 
Among the contributions depending on t — t\, there exist contributions that can increase and decrease 
as the interval t — t\ increases. Having the general expressions at hand, we analyze their dependence on 
the magnitude E of electric field and the time interval t — ti, as well as their behavior at low and high 
temperatures. In the case of a strong electric field, when among the parts depending on t — 1\ the leading 
ones are the contributions due to particle-creation, we are interested in the divergent parts related to a 
large t — t\. In the latter parts, we retain only the contributions leading in t — t\. In particular, we 
analyze the dynamics of mean-energy growth due to particle creation from vacuum and due to the work 
of electric field on real particles. Making a comparison of this dynamics with the energy density of the 
external electric field, we can establish the limits of a correct applicability of the concepts of a constant 
external field. In our calculations, we use the technics and general results outlined below in this section 
with the corresponding citations. 

We assume the initial state to be a state of free in-particles in thermal equilibrium at a given temper- 
ature 9. In such a model, we are going to calculate the mean values of electric current and the EMT of 
the Dirac field. In the general case, the background under consideration is intense, time-dependent and 

2 The physical interpretation in the case of a field that violates the stability of vacuum (for instance, a constant electric 
field) is extremely involved when the field is given by a space-dependent potential which does not disappear asymptotically. 
This happens because the state space of strong-field QFT is quite different from the state space of standard scattering theory, 
in which such a gauge for an external field is appropriate. In the general case of a space-dependent potential, one cannot 
use the standard definition of vacuum in QFT because such a state is well-defined only on a space-like hypersurface. It is 
then generally unclear when a strong electric field is given by a space-dependent potential, what particles and antiparticles 
are, what the vacuum state is, and what is the connection of the definitions being applied with the standard QFT approach. 
In this case, instead of the standard evolution in time, one has to use the rather obscure "evolution in space" . 
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violates the vacuum stability. Such a background must be considered nonperturbatively. In doing this, 
we follow the formulation proposed in [29j . Our calculations can be treated as a one- loop approximation 
within the method of complete QED, developed in [3] [30] and [23l [22] . Namely, this method provides 
a technique of calculating the mean values for systems with unstable vacuum, which is based on the 
generalized Furry picture for systems in strong external backgrounds (see 1291 [3j); this technique may 
be interpreted as a generalization of the Schwinger-Keldysh technique ([31]: see also the review [32]) on 
systems with unstable vacuum. 

It is assumed that in the Heisenberg picture (for notation, see [3J [35]) there exists a set of creation 
and annihilation operators ajj(m), a n (in) of m-particles (electrons), and analogous operators b^in), 
b n (in) of m-antiparticles (positrons), with the corresponding in- vacuum |0, in), and a set of creation 
and annihilation operators ajj(out), a n (out), of out-electrons and similar operators 6^ (out), b n (out) of 
out-positrons, with the corresponding out- vacuum |0, out). By n, we denote a complete set of possible 
quantum numbers. The in- and out-operators obey the canonical anticommutation relations 

[a n {in),a\ n {in)}+ = [a n (out),al n (out)} + = [b n (in), 6j„(m)] + = [b n (out),bl n (out)} + = 5 nm . (4) 

The m-particles (C = + for electrons and C = — f° r positrons) are associated with a complete set 
of solutions {(ip n ( x )} ( we ca ^ ^ an m-set in what follows) of the Dirac equation with the external 
electromagnetic field under consideration with the asymptotics {^ip n (ti, x)} at the initial time-instant 
t%. The latter functions are eigenvectors of the one-particle Dirac Hamiltonian Tt(t) (with external field 
under consideration; see its exact definition in the next section) at t = t%, 

H(tx) c ^(ti,x) = Cs^Mh, x) , H(t) = 7 °([M + 7 (<V - qA (t, x))]), (5) 

where 7 M = (7 ,7) are Dirac's gamma-matrices; q is the particle charge (for an electron q = — e); A 
is the vector potential of the external electromagnetic field (Aq is zero in the case under consideration); 

En^ are the energies of m-particles in a state specified by a complete set of quantum numbers n, and 
s„ > 0. The out-particles are associated with a complete out-set of solutions {^ip n (x)} of the Dirac 
equation with the asymptotics < 'ip n (t2,'X-) at ti being eigenvectors of the one-particle Dirac Hamiltonian 
at t2, namely, 

W(t 2 )^„fe,x)=C4 C)C V'n(i2,x), (6) 

where 4 are the energies of out-particles in a state specified by a complete set of quantum numbers n, 
and 4 ±} > 0. 

The out-set can be decomposed in the m-set as follows: 

<1>(x) = +^(x)G ( + |«) + ^{x)G , (7) 

where the decomposition coefficients C W'j are expressed via the inner products of these sets. These 

coefficients obey unitary conditions that follow from normalization conditions for the solutions. 

The Hamiltonian H (t) of the quantized Dirac field is time-dependent due to the external field. It 
is diagonalized (and has a canonical form) in terms of the first set at the initial time instant, and is 
diagonalized (and has a canonical form) in terms of the second set at the final time instant. For example, 

H (ti) = J2 [4 +) 4 a n (in) + (in) b n (in)] . (8) 

n 

Correspondingly, the initial vacuum is defined by a„(m)|0, in) — b n (in)\0,in) = for every n. 

All the information about the processes of particle-creation, annihilation, and scattering is contained 
in the elementary probability amplitudes 



l + ) 

1 / mn 



w (+|+)„ m = c v 1 < 0, out I a m (out)a'l(in) | 0, in >= G 1 
w(-\-)nm = c v X < 0,out\b m (out)bUin)\ 0, in >= [G" 1 (-1")]^, 
w (01 - +) nm = c- 1 < 0,out \bUin)al(in)\ 0,in >= - [G (+\~) G" 1 (-l")]^ , 
«" (+ - \°)mn = S 1 < 0, out \a m (out)b n (out)\ 0, in >= [G- 1 (+|+) G (+r)] mn , 
c v = (0, out\0, in) , (9) 

where c v is the vacuum-to-vacuum transition amplitude. 

The sets of in and out-operators are related to each other by a linear canonical transformation 
(sometimes called the Bogolyubov transformation). It has been demonstrated that in the general case 
such a relation has the form (see [29] ) 

V (at (out), a{out),b\out), b(out)) = (a f (m), a(iri),tf(iri),b(iri)) , |0, in) = V\0, out) , (10) 
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where the unitary operator V has the forrrJ3 

V = exp {-a\out)w (+ - \0)b^(out)}exp{-b(out)lnw(-\-)b^(out)} 

x exp {a' (out) In w (+|+) a(out)} exp {—b(out)w (0| — \-)a(oui)}, (11) 

and 

c v = (0,out\V\0,out) = exp{-Trlnw(-|-)} . (12) 

Note that the exact formula (fT2|) holds true for an arbitrary external held. In [33], has been demon- 
strated that in the case of a constant electric-like field the Schwinger formula ((T|) leads to the same 
result. 

As our final purpose, we are going to analyze the following mean values: 

(j„ (*)) = Tr [p in j^] , (T„ v (t)) = Tr [ Pm T^} , (13) 

where p in is a density operator of the initial state in the Hcisenberg picture; the operators of current 
density , and energy- momentum tensor T^ u have the form 

U = | Rz), 7m 1>(xj\ , T MV = \ (T— + ) , 

T tT = \ { $(x)n»P» 1>(?)] + [i^(*),7„ } • 

Here, = ic?^ — qA^(x) and VK 2 -) are Dhac field operators in the Heisenberg representation that obey 
the Dirac equation with the external field A^(x). 

In the next section, we calculate and analyze the mean values (|13p for the initial vacuum state. We 
denote them as follows: 

(ip(*))° = <0,*n|i/.|0.*»>. = (0, inlT^lO, in). (14) 

Then, we analyze the more complicated case (|13[) . in which the initial state of the system under con- 
sideration is prepared as an equilibrium state of noninteracting m-particles at temperature 9 with the 
chemical potentials and is characterized by the density operator p in in the Heisenberg picture, 

Pm = Z- 1 exp j/3 ^(Ojv^ - ff (tx) j | , Tr Pm = 1 , (15) 

where Z is a normalization constant, /3 = f? _1 , and wee the operators of in-paricles numbers, 
Af(+) = ^ 4 (in) a n (in) , i\K-) = ^ 6+ (tn) b„ (in) . 

n n 
In our calculations, wc also need the matrix elements 

0m (t)) c = (0,out|jV|0,m)c7 1 , (T MW (t)) c = (0,out|T /u/ |0,in)c7 1 , 
O'm Wont = <0, out|i„|0, out), (T^u (t))° out = (0, outlT^lO, out). (16) 



The matrix elements (jl3|). (|T4|) . and ([16]) can be expressed via relevant singular functions of Dirac's 
fields, as follows: 

*Tr {p in TiP(x)^(x')}=S? n (x,x')+S 6 (x,x') , 
*Tr { Pl J, (x 1 ) i/>(x)} = S+ (x, x') - S e (x, x') , 
zTr {p in ^ (x) V> (x 1 )} = Sr n ( x , x 1 ) + S e (x, x 1 ) , 

S? n (x, x') =0(xo- x' ) S- (x, x') -9(x' - x ) S+ (x, x') , (17) 

where are singular functions at zero temperature, 

s in( x > x> ) = i(0,in\Ttp(x)ij(x')\0,in) , 
Sfa(x,x f ) = i{0,in\tp(x) : ili(x')\0,iri) , 

S+(x,x') = i(0,in$(x')ip(x)\0,in) ■ (18) 



3 Here and elsewhere, we use condenced notations, for example, 

bw (0\-+)a=J2 b n w (0| - +), 



All the singular functions can be expressed in terms of the in-set of solutions as follows: 

71 

S ) Sq ^ ~~~~ Sq (--^5 *C ) 5 

( Xj x ') = £ ±V , n W± ^ (a: ') tv(±) (in) , 



where 



NP (in) = [exp [p ( e «) - pCO) } + 1 



(19) 



(20) 



It is important to emphasize that, in the general case, Sf n is different from the causal Green function, 



S c {x,x') = i(0,ont|TV'(a;)^(a; / )|0 ! m)c- 1 . 



(21) 



From 1[2"T]). it follows that 



5 C (x, a;') = (xq — x' ) S (x, x) — 9 (xq — xq) S + (x, x) 



S- [x 1 x') = iY J + ^ n [*)G 

S + (X,X') =i^2 -4>n ( X 



+ 1 + V>m i X ') , 

/ 71771 

t 



G _ 



i>m i X ') > 

71 

see [25]. Then the difference S p (x,x') = 5f n (ar, a;') — S c (x,x f ) has the form 

^(xy) = ;£ [G( + r)G(_r)- 1 ]^ m+ ^ m (o ; '). 



(22) 



(23) 



This function vanishes in the case of a stable vacuum, since it contains the coefficients G( + |~) related to 
the mean number of created particles. 

Another kind of Green's function that will be used later on is given by 



(24) 



(25) 



Sout( x > x ') = i(0,out\T<ip(x)ip(x')\0,out) . 
It is related to S°(x,x') as follows: 

S c out (x,x') = S c (x,x') + SP(x,x / ), 

S*(x,x') = -iJ2 + i>n( x ) [G( + | + )- 1 G( + r)] nm -^ m ( a; '). 



We can see that there exist relations between the mean values (fl~3|) , (jT4j) and matrix elements (fl"6 
namely, 

0m (*)> = 0m (*))° + (U (t)) 6 , (T^ (*)> = (T»u (*))° + <2> 
& (*)>° = 0m (*)> c + & M) p , (V (f)>° = (t^ (t)y + (t^ (ty, 

(J, (t))° out = & + (J, (t)) P ■ (V (t))° out = (T, v (t)Y + {T, v (t))P . 
These relations involve the quantities 

(j M ( t )) c >™> e = iq tr [^ S °**>°{x,x'j\ \ x=x , , 
{Tftu (t)) c >™> 6 = i tr [A^S c >™> d (x, x')] \ x=x , , 
V = V4 [7 M (Pu + K) + 7 U fa + P'l)] , 



(26) 



P>* = -i— 



qA^x') . 



(27) 



where tr [• • • ] is the trace in the space of 4 x 4 matrices, and the equality x — x' is understood as follows: 

1 



^[■■■( x , x ')}\ x = x > 



lim tr [• • • (x, x )] + lim tr [■ • • (x, x)] 



(\ 



The quantities ([2"T]) are expressed in terms of S c (x, x'), which is the causal Green function (propagator) 
of the Dirac equation with an external field, and in terms of S p,p6 , which are solutions of the same 
equation. The quantities < (t) >°, < (t) >°, < (t) >° ut , and < T^ v (t) >° out are contributions 
to the corresponding mean values at zero temperature and density. The quantities < (t) > e and 
< Tfj, v (t) > e present the contributions due to the existence of the initial thermal distribution. 

The quantities (j M (*))°, (T^ (i))°, (j„ (*)>°„ t > (T^ (t)) a out , (j M (t)) e , and (T^ (t)) e are real-valued by 
construction, due to the properties of the singular functions Sf n , S% ut and S . On the contrary, the 
quantities (j M (t)) c,p,p and {T^ (t)) c ' p ' p are not necessarily real-valued. For the purpose of the following 
consideration, it is useful to rewrite (j M (t))° and (T^ (t))° in the form 

(j M (t)>° = Re(j,(t)y + Re(j,(t)y, 

(T^ (t))° = Re (T^ (t)) c + Re (T^ (t)) p , 

UAt)) out= Re M (*)) C + Re M (*)>?, 

(T, v (t))° out = Re (T, u (t)) c + Re (T, v (t)f . (28) 



3 The initial state as vacuum 
3.1 T-constant field regularization 

In this section, we examine the quantities < (t) >° and < T^„ (t) >° that represent contributions to 
the corresponding mean values at zero temperature. In order to calculate such quantities, the T-constant 
field regularization is necessary. In particular, we are going to study the leading contributions to the 
divergent parts of the quantities at T — > oo. To make the consideration complete, and to provide the 
reader with some formulas necessary for the further consideration, we start by reproducing some of the 
results on the mean numbers of created particle that were obtained in [28] . 
We chose the T-constant field potentials as follows, Aff = Af = A® = 0: 

(Eh, te (-oo, ii), 
Af(t) = { Et, te [h,hl 

{ Et 2 , te {h,+oo), 

where t% = —h = T/2. The corresponding electric field E(t) is given by 

r 0, te (-oo, h), 
E(t) = { E, te [h,h], (29) 

{ 0, te (t 2 ,+oo), 

First, let us suppose that the magnetic field is absent. 

For the purpose of T-regularization, it is sufficient to choose, for {AAA} and A^AA} , defined 
above, the following orthonormalized sets of solutions of the Dirac equation with a constant electric field: 

±Vp,r(aO = il P + M )±4> P ,±lA X ) > ^ P A X ) = ^ P + M )^P,Tl,rW . 

±^p,,,rW = ±^p, s (i)exp{ipx}u s>r , ± <t> p , s A x ) = ± ^ P , s Wexp{ipx}w SjI ., 

+<f> Pta (t) = CD V _ i+. (±(1 - i)0 , i0 p , a (i) = CD -»-^ (±(! + 00 ,s = ±l, (30) 

where p is the momentum and r = ±1 is the spin projection; D v [z) is the Weber parabolic cylinder 
(WPC) function [34], and 

iX M 2 + p 2 ± j 2 qEt - p 3 

u =-^> A = — I F , A Pi. = (p ,P ,0), C = = — 7f=FT> 

C - (2tt)- 3/2 (2 l^l)" 172 exp (-ttA/8) , 

w Si r are constant orthonormal spinors, v\ r v s y — 5 r y, subject to the supplementary condition (1 ± 
7°7 3 )i> T i ir = 0. Using an asymptotic expansion of the WPC-function [53], 

A,(*)=^exp(-£) (E H ^l~ 2) t )n + °(\*\- 2{N+1) )) » i ^ < W4, (3D 
we can obtain the asymptotics of particle energies as T — > oo, 

4° = I^+P»I, ^ = l^-ftl- (32) 
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One can see that the matrices G ( ^ J ([7|) are diagonal: 

G(cl C ') ,=^,,'<5(P-P')9(cl 

V / p,r,p ,r ^ 



(33) 



The differential mean numbers of electrons (equal to the corresponding differential mean number of 
pairs) with a given momentum p and spin projections r created from vacuum are 



51-1 



(34) 



<5p,p'- Note that H p r is even with 



^p,r = (0 7 in\a^ r (out)a Ptr (out)\0,in) 

were the standard volume regularization is used, so that <5(p — p') 
respect to the sign of qE. 

If the time T is sufficiently large, T >> To = (1 + X)/y/\qE\, the differential mean numbers H PiI . have 
the form 

3^ 



1 + 



1+A 



O (1) , -K p < £ < K p , 

3 N 



(35) 



O 



1 + A 



where 



\ P3 \-\qE\T/2 



and K p is a sufficiently large arbitrary constant, K p >> 1 + A; see [28]. In the limit T — > oo, the 
differential mean numbers have a simple form: 



H P) , = e- 7rA , 



(36) 



which is identical with the one obtained for a constant electric field by Nikishov [33J . One can see that the 
stabilization of the differential mean numbers to the asymptotic form (|36[) for finite longitudinal momenta 
is reached at T >> To- The characteristic time To is called the stabilization time. 

In order to study the effects of switching the electric field on and off, at T >> To, one can study 
another example of a quasiconstant electric field: 



E(t) = £JcosbT 2 ( - 
1 a 



(37) 



This field switches on and off adiabatically as t — > ±oo and is quasi-constant at finite times. The 
differential mean numbers of particles created by such a field have been found in [35] . As shown in [28j , 
the differential mean numbers in the field (|37| have the asymptotic form (|36p for sufficiently large a, 
a » cto = (1 + VA)/VW and for \qE\ 

a 3* \p3\- Thus, ao can be interpreted as the stabilization 
time for such a field. At the same time, the latter implies that the effects of switching on and off are not 
essential at large times and finite longitudinal momenta for both fields. Extrapolating this conclusion, 
one may suppose that in any electric field being quasiconstant ~ E at least for a time period T » To 
and switching on and off outside this period particle-creation effects have no dependence on the details 
of switching on and off. Therefore, calculations in a T-constant field are representative for a large class 
of quasiconstant electric fields. 



3.2 Contributions due to particle-creation 

Let us apply the above results to a calculation of the leading terms in Re(j /i (t)) p '^ and Re(T M „(i)) p 'P 
from ([28]) for sufficiently large times t, more exactly, for large values of t — 1\ =t + T/2, and of T. Being 
written in a dimensionless form these two conditions have the form 

^\q~E\(t + T/2)^>l + M 2 /\qE\ , t<T/2, 

yf\qE\T ^> I + M 2 / \qE\ , t>T/2. (38) 
We note that the dimensionless parameter 



VW\T = 




that has emerged here will play an important role in the further considerations, and will enter all the 
consistency conditions that are obtained in this article. 



x 



We can see that ([55)1 is a stabilization condition when the leading terms do not depend on the details of 
switching on and off. As follows from ([77)1 , we need the singular functions S p (|2"3")l and S p at x w %' in 
such an approximation, which provides the required contribution for Ke(j^(t)) p ' p and He(T tlv (t)) p ' p when 

(|38|) is valid. Time-dependence arises due to integration overp3 under the condition — y/\qE\-j < f < —if 
with the subsidiary condition that |£(t)| must be sufficiently large, |£(t)| > K, where K is subject to the 
condition K ^> 1 + to 2 / | Q-E j . The distribution H p r plays the role of a cut-off factor in the integral over 
p_L, the T-dependent contribution of S p ' p being thus convergent. 

The range of integration over momenta in S p that determines the leading contributions can be defined 

as 

D .i |p-l| < VWllVWlit + TW-K} 1 ' 



-T/2 + If/VM < p 3 /g£ < t - K/y/\qE\ 

Using to express the solutions ±"0n y i a ± '0ni an d taking into account the asymptotic expansion 
(J2J), we can calculate the leading contributions to the function S p (x,x') (j2"5j) . 



/ dp V Hp . r [+V' Pir (a;) + 7/; p . J ,(a;') - 4> Ptr (x) 4> p , r (x')] 

J D P=± i 



Taking summation over r, we represent (|39|l as follows: 

S p (x,x') = (jP + M)A p (x,x'), 

A p (x, x') = i [ N p , exp{*p (x - x')} [+0 P , O (i) + ^, (*') +" ^oW^oCO] d P • 
Then, integrating over in (|4"0")) , we obtain 

, h\\[x\i,x»)h±[x.±,x ± )d — , 

-T/2+K/ V / j^E| V 7 V^/ 



J-'_ 



where 



^rJ^exp{i[^) 2 -^) 2 ] + iP3(,3-4) 



/i(x ± ,x' ± ) = (2vr)- 2 |g£;|exp - 



(x x -x' x ) 2 | g £| 



1^1 



4tt 



and the notation xtf = 0, 0, x 3 ) and xj_ = (x 1 , x 2 , 0) has been used. 

In the same manner, we calculate the leading contributions to the function S p (x,x') 



S p (x,x') 



( 7 P + M)A p >,x'), 



A p (x, x') = i H p ,o exp{zp (x - x')} [ + ^ (t) + <f>* pfi (t / ) tfp, (tMp,o(*')] d P • 



D 



The range of integration D reads 



|p±| < VkE\[VW\(t + T/2)-K 
t + K/^\qE\ < P3/qE < T/2 - K/y/ 



1/2 



qE\ 



Integrating over p± in (|43|) . we find 



A p (x,x') = i 



T/2-K/J\qE\ 



t+K/J\qE\ 



h\\ {x\\,x'^j h± (x±,x' ± ) d f^J , 



(39) 



(40) 



(41) 



(42) 



(43) 



(44) 



The following results will be outlined in the presence of a constant magnetic field, B ^ 0. For such 
a field, we select the nonzero potential as Af = Bx 2 . Now, the complete set of quantum numbers that 

describes particles is (pi, Ub,P3, t), hb = 0, 1, . . ., A = (M 2 + \qB\(2riB + 1 — r)) \qE\ . The space-time 



Q 



dependent part of the function <j)p,s,r( x )i see (|50|) . is modified as follows: 



(27T) 



" 3/2 exp{ipx} -> (2tt) 1/2 exp{-ip 3 a; 3 }(/) pi! „ B!j ,(x ± ) , 



.(x±) = 



2 n B +l 7r 3/2 n I 



exp < —ipix 



\QB\ 



Pi 
qB 



(45) 



Here, 7i nB ( x ) are the Hermite polynomials. Then, integration over pj_ in the above formulas must be 
replaced by integration over pi from — oo to +00 with summation over the integer quantum numbers ns 
in the interval 

n B < [VW\ (t + T/2) - K\ \E/2B\ . 

After this, one can see that formula (|41[) needs only one modification. Namely, the new (in the presence 
of magnetic field) value of h± (xj_,Xj_) reads 



h± (xj_,Xj_) = 



qB exp(7rE 3 B/£;) 
47rsinh(7rS/E) 



■ exp 



{-^-(xx-xl) 2 fcoth ( ^)} 



(46) 



Using and JSj), one can represent (j M (<)) p >? and (T M „ (<)) in ((27]) as follows: 

(j„ (*)> p ' p ~ - tr ^PA^Car, *')] | x=x , , 



<V (t)>™ 



tr ^V 7 PA^>,z')]| x=a , 



Then, using (gTJ, gH), JUJ), and (gBJ), we find 

OmM) p ' p = -^P^trA^^ar')!^, 
<T w (i)>^ = -iPfcA p > p (x,x>)\ x=x , . 



(47) 



The off-diagonal matrix elements of (T , At „) p and (T^ u ) p are all equal to zero. 

Taking derivatives, calculating traces, and integrating over p 3l we obtain the leading contributions at 
large T. First of all, 

(3, (t)) p = -2%q sgn (qE) (1/2 + t/T) 
Om (*)>* = -2^9 sgn (1/2 - i/T) " cr , 
(Too (i)) p = (T 33 (t)) p = \qE\ T (1/2 + t/T) 2 n cr , 
(Too (t)) p = (T 33 (t)) p = \qE\ T (1/2 - t/T) 2 n c \ 



where 



Secondly, 



4tt 2 



EBT coth(ir B/E) 



exp — 7T 



M 2 

W\ 



o 



\VW\t / 



(48) 
(49) 



(Tu (t))" = (T 22 (<)) p , (Tn (<)) p = <T 22 (t)>* 
(T n (t)>* = n f ln 

(T n (t))P = n 



y/\qE\(T/2 + t)\ +0(]nK), y/\qE\ (T/2 + t) > K 
0(\nK) , ^\q~E\(T/2 + t)<K 

-ln 



VPf(r/2-t)] +OQnK), y/]qE\(T/2-t)>K 
-0(\nK) , ^(T/2-i) <K 



where 



(qBf 



A-k 2 sinli (ttB/E) 



exp — 7T 



Af 2 



(50) 



(51) 



We can see that all the leading contributions are real-valued. Note that the current density (j^ (t)) p 
and the component of EMT (T^ u (t)) p are zero for t < t\ and change with time until i 2 , (j 3 (t)) p being 
linear, (T 00 (t)) p = (T 33 (t)) p quadratic, and (Tn (t)) p = (T 22 (t)) p logarithmic, respectively. As will be 



in 



is the total number-density of pairs created by the T-constant electric 
field. 

After switching on of the electric field (at t > t 2 ), all the mean values (|48|) . (|50|) are constant and 
retain their values at t 2 . 

We should mention that formulas (|48|) . (|50|) have been obtained for the first time. 
At large times t — t% = t + T/2 ^ K/ y/\qE\, it is useful to examine the quantities 



fl (*) = (3, (*))° - <j„ , TZ (t) = {T, v (t))° (T„„ (i)>°, 



(52) 



which can be interpreted as those corresponding to oui-particles at any large t. At the final time instant 
t 2 , they coincide with the physical quantities (j^ (t2)) p and (T^ (t 2 )) p ■ Then it follows from (f2"6"]l and 
ffl, that 



and 



(«) = (?„ W) P ~ & (t)> ? = -2^«sgn( 9 £0 (2i/T)n c 
T£ (t) = (T lw (t)Y-(T, u (t)) p , 



(T 00 (t)r r = (T 33 (t)) cr = 2\qE\tn cr , (T u (t)) cr = (T 22 (t)) cr , 



(53) 
(54) 



(T U (t)) cr = 




\ q E\([T/2f -t 2 )] +0(lnK), ^\q~E\ (T/2 - t) > K 
(T/2 + 1)]+0 (In K), ^/\q~E\(T/2-t) < K 



At t > t 2 , the production of pairs terminates, and the quantities ([53"]) and (JHH) maintain their values 
(|53|) and ([54]) at i = t 2) and present the current density and EMT of created particles. For example, 
from j c J~ in (|53[) we can see that n cr is the total number-density of pairs created during the entire time 
of action of the electric field. 

We emphasize that expressions ([4*5]). (IBTfl) for {j^(t)) p and (T fiu (t)) p represent contributions to (j M (i))° 
and (T^(t)} , respectively, at any time instant t if the stabilization condition (f3"5| is valid. On the 
contrary, expressions (f5"3"]) , ([54"]) are valid only for time instants i that are sufficiently close to t 2 , 

i: (T-2i)/T< 1, 

when the interpretation in terms of final particles already makes sense. 



3.3 Vacuum polarization contributions 

We are now going to calculate the real- valued parts of the quantities (j M (t)) c and (t)) c defined in 
(j2~7| . which do not have any T-divergences. In these expressions, we can use the causal Green function 
S c (x,x') (|22[) in the constant field, since the above quantities do not need any regularization. 

For such a Green function, the so-called Fock-Schwinger proper time representation holds true: 



(55) 



S c (x,x') = (jP + M)A c (x,x'), A c (x,x')= f(x,x',s)ds; 

Jo 

see [33] and [36], where the Fock-Schwinger kernel f(x,x',s) reads 



f(x,x',s) 
fW(x,x',s) 



exp 



q 2 EB exp (iqA) 
(4ir) 2 sinh(qEs) sm(qBs) 



exp 



-iM 2 s - i 



(x — x 1 ) qF coth(gFs) (x — x') 



(56) 



see [5]|3Z]. Singularities of the kernel are situated at the origin and on the imaginary axis, \qE\ s = —iim, 
n = 1,2,3,... . 

It should be noted that in this expression the only term 



A = A||+A_l, A ll = {x +x' Q )(x 3 ~x' 3 )E/2, A x = - f A^dx^ , 

J x' 



(57) 



is potential-dependent. Here, is the potential of magnetic field F^, and the integral is taken along 
the line connecting the points x and x' . 

Using this representation, we calculate Re(j M (t)) c and Re(T M „ (i)) c - It is easy to see that (j M (t)) c = 0, 
as should be expected due to translational symmetry, and (T^(t)) c — , ji 7^ v . 



1 1 



Calculating the diagonal terms (T^^it)} , we discover that it can be derived from the Heisenberg-Euler 
Lagrangian £, 



1 f°° 
£ =—tr / s f(x,x,s)ds . 

2 Jo 



Subtracting the zero external field contribution from £ and performing the standard renormalizations 
of £— £\ F=0 , that leave qF^ invariant, i.e., the standard renormalizations of the coupling constant q 2 
and potentials A^, we obtain a finite expression: 



ds exp (— «M 2 s) 
8^ 



.,2 



I a 

q 2 EB coth (qEs) cot (qBs) ^ — — (E 2 - B 2 

s z 3 



(58) 



see [7|. Making the same renormalization for (2^(4))°, we can see that there holds the relation 

{Too{t))ren = ~ (^Wren = E dE^ ~ ^ ^ ren ' 

(Tn(t)) c ren = (T 22 (t)y ren = £ (t) ren - B 9C{ ^ , (59) 

where £ (t) ren = £ re n\ E ^E(t) > and E ( t ) is defined b}0 fl25J|. 

Finally, the vacuum mean values of the current density and EMT have the form 

(j M (t))° = Re^t))*, <T M , y (t))° e?1 = Re(T^(t)) c ren + Re^t))*, (60) 

where Re(j M (t)) p and Re(T Ml/ (i)) p are presented by ijlSjl and ([50)) . We can see that the T-dependent 
contributions Re(j fl (t)) p and Ke(T /1 „(t)) p arise due to vacuum instability; they are global physical quan- 
tities and have the factor exp {— ttM 2 / \qE\}. This factor is exponentially small for a weak electric field, 
M 2 1 \qE\ ^> 1, and these quantities can actually be observed as soon as the external field strength 
approaches the critical value E c . On the other hand, the term Re(T A1!/ (t))J; en is t-dependent and T- 
independent; it is therefore local and does exist in arbitrary electric fields. When the T-constant electric 
field switches off, the local contribution from the electric field to Re(T^(t))^ en vanishes, but the global 
one, given by Re(T Aljy (<2)) p , does remain. Thus, in the general case, both kinds of contributions are im- 
portant. We stress that in order to make (2)«,(i)}° en finite the only term {T^, v {t)) c has to be regularized 
and renormalized due to the standard ultraviolet divergences, which is consistent with the fact that the 
ultraviolet divergences have a local nature. 

Using expression (|60p . we can find a condition which imposes restrictions on the strength of classical 
constant electric field. Note that for a strong electric field (B — 0), M 2 / \qE\ <C 1, and for a large T, 
there is a well-known asymptotic expression for the vacuum energy density Re(Too(t)}^ en , 

Re(T Q0 (t)y ren = -^E 2 Jj^. (61) 

It does not depend on T. The energy density of a classical electric field is E 2 /8ir. Then (see [55]) one 
makes the conclusion that the notion of a strong constant electric field is physically meaningful when 
~TT|) is much less than E 2 /8w, which implied 



However, if T is large, one also has to take into account the T-dcpcndcnt term Re(Too(i)) p in (f60|) . 
At t = t% , and for an arbitrary magnitude of the electric field, we have 

Re (T 00 (t 2 )) p = ?A \qE\ T 2 exp ■ (63) 



4 In the absence of electric field, the Green function S c (x,x'), Sf n (x, x'), and S^ ut (x, x') coincide. Therefore, in case 
t<h we have {T^t)) =. (T M „(t))° \ E=Q = <T M „ (i)} c | B=0 . For t > t 2 , relation JUJ implies 

(T^(t))° = (T^(t))° out +T^(t) , 

where T^ v (t) is given by J53t.ll54t and {Tfj. v (t))° ut =. (T MI/ (t))° \ E _ Q = {T^v(t)) c \ E=0 , i.e., this is a vacuum polarization 
contribution. 

5 The same inequality has been obtained by Ritus as a condition to avoid the breakdown of the loop expansion in a 
strong field; see, e.g., 1391 . 
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One can neglect the back-reaction of created pairs in electric field only in case ([55)) is far less than E 2 /8ir, 
which implies 



or, in Gaussian units, 



|^|^«^ex P ^), (64) 



C , , 7T 2 C fi / C 3 M 2 

-\qE\T<^ — — exp 7T 



1 2q 2 h \<l E \ 

The restriction (|64|) is completely new, it is much more restrictive than (|62|) . ft imposes restrictions 
both on the maximal strength of electric field and its duration due to the back-reaction of created particles. 

4 Initial state as thermal equilibrium 

We are now going to calculate the quantities < (t) > and < T M „ (t) > (fl3|) with the density operator 
(fi"5"|) . by using relations (f26|) . In fact, the contributions < (i) >° and < T M „ (i) >° to these quantities 
have been calculated in the previous section. Here, we study the remaining temperature-dependent 
contributions (j M (t)} and (T M „ (t)) e to these quantities. 

4.1 Contributions due to the particle creation 

Let us examine the function S e (x, x') in (fH)]) . The form of the functions Sj in (flU)) for t < t\ , and i' < £j 
is known, they are functions for zero electric field. Therefore, we have to study the case t > t\, t' > t\. 
First of all, we separate the contributions from particle creation in the same manner as it has been done 
for the function Sf n (x,x') in Sec. 2; for details, see Subsec. 1 of the Appendix. 

Using (|26|) . ([60]) . and (|136p (see the Appendix), we present the mean values (j^ (t)) and (T /tv (t)) ren : 

(3, (*)) = Re(v (*))§ + J£ (t) , 

J»(t)=ReOV(*)> p + BfiO'/« (*)>?, 

(T^(t)) reTi = Re(V(t))£ e „ + Re(T„„ (t))g + (t) , 

7*, (i) = Re^f))* + Re(T M , (t))f . (65) 

Now, {Tui/(t)} is replaced by the renormalized (with respect to ultraviolet divergences) quantity {T^ u {t)) ren 
by inserting ^r^ v {t))^ en in the r.h.s. instead of (T^t))". The terms J? (t) and t^ v (t) with the upper 
script p represent the contributions from particle-creation, the terms with the upperscript c represent the 
remaining contributions. One ought to recall that the vacuum polarization contribution {T^vit))^^ has 
been already calculated and has the form (|59|) . In the absence of electric field (for t < t± and t > fa), 
the quantity (|65|) does not depend on time, in particular, JP (t) — t\ v (t) = for t < t%. As has been 
already mentioned, the interpretation in terms of out-particles makes sense only for t > t, where t is 
sufficiently close to t 2 , (T — 2t) JT <C 1. At the same time, for t > i the quantities J? (t) and t p (t) are 
the current density and EMT of out-particles created by the electric field. The quantities Re(j M (t))g and 
Re(T„ I/ (t))g describe the corresponding contributions related to the existence of real in-particles at the 
initial state. 

The quantities Re^ {t)) p g and Re(T M „ {t)) p e present the contributions from particle-creation at finite 
temperature. They are T-divergent. The leading contributions to these quantities at large T can be 
found by analogy with the case of zero temperature, which has been investigated above. We examine the 
leading contributions to Sq which are given by the expression 

s p 6C (x,x') = i f dpY, (*»)*W ( + %A X ) + %A X ') - ~%A x y%A x ')) • ( 66 ) 

Jd r=±l 

The integral J D dp is defined in (|3D[) : the differential numbers H p r are defined in (|3"6")l . and 

ATjO (j„) = {exp [(3 (e« - //<«)] + (67) 

where £p in ([67]) are quasi-energies at t±; see ((32)) . 

In order to obtain a generalization of (|6"6"|) to the presence of a constant magnetic field, B ^ 0, we 
have to follow the way described in Subsec. 3.2. Then, taking summation and integration, we obtain the 



13 



leading contribution to the function Sg ^ (x,x') in the case under consideration: 
S p e c (x,x') = ( 7 P + M)Al c (x,x') , 

r t-K/J\qE\ 



AlAx,x') = 



T/2+K/J\qE\ 



h± (x±,x 1 )d 



P3_ 

qE 



(68) 



The functions h\\ (x\\,x'^j and h± (xj_>Xj_) are given by expressions (|^2")) and (|36")> . respectively. 
Therefore, we obtain the leading contributions to (j^ (t)) p e and (T^ (t)) p 9 , namely, 



* ^ ^ \x—x' 



(T„»(t)) p e = 0, il±v; (T MM (t)) e 



-^P^trA^z') 



(69) 



In order to examine the temperature-dependent contributions to J p (t) and t p (t) in (|65[) , we have to 



investigate the quantities ([69]) . At low temperatures, ft (ep" 1 — /Li(0 



3> 1, and assuming £jP > \u 



i.(0| 



one can see that the contributions from Re(^ {t)) p g and Re(T My (t)) p g are very small in comparison with 
the vacuum contributions Re(j M (i)) p and Re(T /J „(i)) p . Therefore, in the case under consideration the 
leading contributions to J p (t) and r p „ (t) are given by 



J p (t) = Re(j,(t)) p , r p (t) = Re(T^(t)) p . 



(70) 



For £p < 1^(01 with a sufficiently large \fi^\ 3> y/\qE\K, contributions from Ke(j^(t))g and 
Re(T M „ are comparable with contributions from Re(j /J (i)) p and Re(T M „(t)) p , respectively. For ex- 
ample, at |/i(0| > (f + T/2) we have 



Then 



Ml* (t))g = -2Re<j p (t))*, Re(r p „ (i))f = -2Re(T^(t)) p . 



J? (t) = - Re(^)) p , r p (t) = - Re(T^(i)) p 



(71) 

The negative sign in these expressions implies that the total number-density of particles decreases due to 
electron-positron annihilation in electric field, and, as a consequence, (j^ (t)} and (T^ w (t)) re „ in (|65[) arc 
less than Re(j M and (7/^,(4))^ + Re(T M1/ (i))g, respectively. This fact is in agreement with the fact 
that, due to particle creation, there holds the relation 

N$ (out) = (1 - H m ) Njp (in) +H m [l - N^(in) 

between the initial Nm\in) and final N^p (out) differential mean numbers; see [22]. Thus, the differential 
mean numbers of created particles are given by the difference ANm = Nm^ (out) — Nm (in) . As a result, 
AN^ is negative in case N£ ''(in) + N^\in) > 1. 

At high temperatures, ft [\qE\ (t + T/2) - fi^] < 1, we find 



\ [\qE\ (t + T/2) - M (+) - } + (if) 



Re(j„(t)> p : 



T?l(*)=^ 2 (*) = § 



i (i + T/2) - \ (,/+> + + (X) 



Re(T 00 (t)) p , 



| g £|( f + r/2)--( M (+)+ / x ( - ) 



1+ ^\q~E\0 



In if 



In 



^^(i + T/2) 



Re(T n (i)) p . 



(72) 



We can see that all the components of J p (t) and t?„ (i) at high temperatures are far less than the 
corresponding components of the vacuum contributions (j^(t)) p and (T^ v (t)) p , which is natural for Fermi 
particles. 
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The resulting expressions for the mean energy density of created particles make it possible to conclude 
that the restriction (l6"4")) for the external constant electric field, established for vacuum, is also valid at the 
initial low-temperature state. In case the temperature of the initial state is sufficiently high in comparison 
with the mean kinetic energy of created particles, (3 \qE\ T <C 1, this restriction changes considerably. In 
this case, the mean energy density of created particles after the electric field turns off, 

^oo (*2) = \p\qE\TRc(T m {t 2 )Y , (73) 

in accordance with ([72]) . is much smaller than the corresponding vacuum expression Re(T o(£2)) p in ([63]) . 
Repeating the arguments of Subsec. 3.3, we find that one can neglect the back-reaction of created pairs 
in electric field only in case Tq (£2) is far less than E 2 /8tt, which implies 

^' T3 «f°4'm)- (74) 

This is much less restrictive than 



4.2 Conventional contribution 

We are now going to calculate the temperature-dependent contributions to the mean current density 
and EMT, which are labelled by the upperscript c. We refer to these contributions as conventional 
contributions. Such contributions are not related to either particle creation or vacuum polarization. 
Thus, it is natural to say that they are determined only by the quantum statistical behavior of the initial 
thermal gas of charged particles in the constant external field under consideration. 

4.2.1 Proper-time representation 

Let us represent (j^ {t)) c g and (T^ v (t))g in (|65|) as integrals over the proper time. To this end, we relate 
S$£ in (|135p (see the Appendix) to the proper-time representation {55)) of S c ; for details, see Subsec. 2 
of the Appendix. For simplicity, we assume that the chemical potentials are not very large, < M. 

As always, we consider large T, namely, 



\f\qE\T > 1 + M 2 / \qE\ . 



(75) 



Using (I143p . Q145|) and (|147p . as well as the explicit form (|151[) for g(x,x' , s,t), and calculating all 
the derivatives and traces, we find the following expressions for Re(j M (t))% and Re(T MI , (£))#: 

Re<i,(£) >g=£ £(-!)<+ Re jf, 
C=± 1=1 



J, W = 



h (s, u) 



du h (s, u) Y (s, u) jfj, (s, u) , (s, u) = iSlq 2 E (t + T/2) 



q 2 EB cot [qB (s - ir)} 



87r 2 V7r u 1/2 sinh (qEs) 
Y (s, u) = exp 



1/2 



e -u/4 e -M'(is+r) 



(qEf{t + T/Zf — 



(01) In, 



and 



(76) 



Re(^(£))g = ^^(-l) i+1 e ^ (0 ReTW, 
C=± l=i 



r (0 



ds 



du h (s, u) Y (s, u) t^v (s, u) 



where it is only the diagonal values £ M „ (s,u) that are nonzero: 

t 0Q (s,u) = ^ t kk (s,u) + M 2 cosh(qEs) , t%x (s, u) = £22 (s, u) 

i {qEf{t + T/2f 



-iqB 



fc=l,2,3 



sin [2qB (s — ir)] 



(s,u) 



2 (a — ir) cosh (qEs) 



cosh (qEs) 



(77) 



(78) 
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Note that the integrals over s and u in the representations (|75)l and (|77|) are finite. Such representations 
hold true for any E and temperature for t £ [£1,^2] within the restriction (|75p on the interval T. 

Representations (1761) and (|77|) could be transformed to other forms, being more convenient for regions 

near s = in the integrands; for details, see Subsec. 3 of the Appendix. The result is that Re J^' and 

Re Tjiv in ([76]) and (f77|) can be presented as 



where 



and 



Re J« = J° + Re J* , Re T« = f « + Re T* , 



TO 



rpa 



Ttl\2qE\-Q 



ds' 



ir/\2qE\~0 



ds' 



du'h{s',u')Y(s',u')~i li {s',u') , 

JC 

<fet' h[s',u')Y[s',u') V (s',it') 



/i (s', u') = —ih (i (cq — s') , v! + uq) 
1 q 2 EB coth [qB (s' + St)} 
~ 8n\/n u' 1 / 2 S in (g£ s ') 6Xp 

F (s', u') = y (i (cq — s') , u' + uq) = exp 



M [s' + <5r) 



(g£) 2 (i + T/2) 2 tv! 



(u' + uo) 1 + 0-Et)* 



St = t- 1 arctan [qEr) , r = -7 

7 / ' 



u + u 



(79) 



(80) 



(81) 



j fl [s',u') = 

%{s',u') = 

i n (*', u') = 

^3 (*',«') = 



(* ( c o - s ') ! u ' + u o) , ?/n/ (s', u') = ^ (i (c - s') , u' + uo) , 
3 q {qEf (t + T/2) r [1 + g£r tan (gSa')l 



t 22 (a 7 , «') 



tan (g^s') 1 + (g£r) 
" sinh [2g5 (s' + <5r)] ' 



qE 



{qEf [t + T/2) 2 t 2 [1 + qEr tan [qEs')] 



2 sin (g.Es') 1 + (qEr) 



1/2 



sin (g.Es') tan [qEs') 1 + (gTr)' 



3/2 



— , 1 i . ~ / ; ;s . r? / t-, /s 1 + ff-ET tan [qEs') 
t 00 [s',u') = 22 t kk [s',u') + M 2 cos [qEs') 2 y -± TJ j± 



fe=l,2,3 



1 + [qEr )' 



(82) 



Note that r < (g_E) 1 tan (gi£s') in and . The expressions for and are given by (|153|) and 
(|154[) , respectively. The form (|T9[) is convenient for the study of both low-temperature approximations 
to a weak electric field and high-temperature approximations to an arbitrary electric field. 

In a weak electric field, \qE\ /M 2 <C 1, the quantities Re and ReT* are exponentially small, so 



that in order to find the expansions of Re jjp and ReT^'J in the powers of \qE\ /M 2 ,we use the relations 



(0 w 



where the exponentially small contributions have been neglected. 



(83) 



4.2.2 Zero electric field 

When the electric field is zero, the mean current density is zero as well. Considering the limit E — > 0, one 
can see that the integral (|80|) over s' for T^ v is taken from zero to infinity. In addition, the dependence 
on the variable u' in the kernel of the matrix elements T" factors out as the multiplier e~ u / 4 /y/u' '. 



1fi 



Calculating the integral over vl by using (|139[) , in case E = 0, we find 



TL=TL= I dsp (s)tl,(s) 



qB coth (qBs) 
Po ( s ) = 4^ ex P 

*?i 00 = 4 00 = 



4s 



2s 



sinh (2qBs 



fc=l,2,3 



The quantities T° are real-valued. Integrating by parts, one can see that Xq has the form 



(84) 
where 



tfn 5 =- — 



2s 4s 2 

In the case under consideration, the vacuum is stable, while the renormalized EMT (|65p is real- valued, 
time- independent, and is given only by the conventional terms 

CV(«)> ren| B=0 — 

enl_B=0 



C=± 1=1 



(85) 



where (^(^enl^o is determined by ^ at E = 0. 

The system under consideration is a relativistic fermionic gas in an external magnetic field, B, im- 
mersed in a sufficiently large quantization volume V. This system is described by the partition function 



Z = Trexp 1/3 I (jl^N^ - H 

that follows from (fl~5|) . H = H (ti), and by the thermodynamic potential fi = ~9\nZ. At the state of 
equilibrium, the temperature- and field-dependent effective Lagrangian C ren is related to O as follows: 



C 6 ren = —Sl/V. 



(86) 



The mean values of energy density (H)/V and particle density in the state (| 15[) can be expressed in 



terms of C ren as 



V 

(C) _ UJ ~ren 



a/? 



c=± 



(87) 
(88) 



Note that the energy density has been already calculated as a component of (|85p . namely, 

(H Q )/V= (T 00 (t)) r enU = o ■ 

Thus, we find that the effective Lagrangian has the form 

f-'ren = ^rm \ E=Q 4~ A£g , 



-M 2 s 



4s 



g£> coth(gi?s) , 



(89) 



where £ren|£ = o is the renormalized Heisenberg-Euler effective Lagrangian (|5"8")l at _E = 0. 

It turns out that the space components of {Tik(t)) ren \ E=0 (|85[) can be derived from the effective 
Lagrangian C e ren as follows: 



(^ll(*))renU=0 _ ( T 22{t))ren\ E=0 — £■ 
(T33(t)} rerl \ E _ = C- ren - 



_ r> ren 

~ ~777T 



(90) 



1 7 



It should be noted that the finite temperature and density effects are contained only in ACg. In 
particular, the thermodynamic potential Qf ree of perfect free Fermi gas at B — is related to ACg\ B=0 
in the form f2/ ree = — ACe\ B=Q V . Expressing the proper-time integral for ACg\ B=0 in terms of the 
modified Bessel function of the third kind (Macdonald's function) K v (z) [34], and using formula 3.471.9 
in [40], 

x u-i e -a/x-b Xdx = 2 ^ay/ 2 R ^ ^ ^ 

we find 

Thus, we obtain for the quantity f2/ ree the expression well-known from textbooks. 

Considering the usual case — — /i^ - - 1 = /x, we can sec that ([89)) is identical with the results 

obtained earlier in the case E — 0: see Eq. (52) in the paper [15] . where the imaginary-time formalism 
was used, and Eq. (5.20) in [16] . where the real-time formalism was used (we disregard the factor 1/2 in 
Eq. (5.20) of [12] as an misprint). In the case /i = 0, Eq. ([59")) is in agreement with the result obtained 
in [II]. 

It should be noted that our analysis deals primarily with the effects of electric field at finite initial 
temperatures, and therefore we do not examine in detail the particular case of magnetic field (at E = 0). 
Note, however, that the leading magnetic field-dependent terms follow immediately from ([59)1 . We present 
these terms for comparison with a nonzero electric field: in the weak-field limit, when either \qB\ /M 2 -C 1 
or \qB\ 1 < 1, we have 



d = AC 6 AC 6 \ B= , = <gf £ £ (-1)< +1 K (M(3l) . (93) 
In the strong-field limit, \qB\ 3> M 2 ,[3~ 2 , we have 

The latter quantity is small in comparison with the corresponding vacuum polarization term, which is 
given by the following expression: 

f | r , {qBf ^ \qB\ 

L -ren\ E = a J -ren\ E=0iB=Q — ^f2 ' 

In the weak-field limit, we can see that the leading field-dependent thermal contributions to EMT 
have the form 



r B — r p dA _ r B 

r 22 °-qq - ->~1 i T 33 



•a - - £ >-^ + g" <0 ^- <94 > 

Note that Fermi systems with an external magnetic field have been under intensive study, and various 
asymptotic expansions of the effective Lagrangian at finite temperatures have been made: for a review, 
see [TSJ |TB] ; a review on finite-temperature free relativistic systems can be found, e.g., in [521 |4"T1 W2\ . 



4.2.3 Weak electric field and low temperatures 

In a weak electric field, \qE\ /M 2 <C 1, the quantities Re Jjj/ and ReT^'J in ([76]) and ([77)) are given by 
and f* v d80))-([82)); see ([83]). The dependence of expressions ([80])-([82)) on the parameters E and (3 is 
described by qEr, which is a function of two variables: s' and u' . Approximations of this function and 
of the corresponding integral over v! at low temperatures and at high temperatures are quite different. 
Therefore, unlike the case E = 0, these limiting cases must be studied separately. In this subsection, we 
shall examine the case of low temperatures, M/3 ^> 1. 

Contributions to integrals ([80)) . given by the region M 2 s' 3> 1 of large s', are exponentially small. 
Therefore, we limit the range of integration over s' from above by the value x/M 2 , so that 1 C x ^ 
M 2 1 \qE\. At low temperatures, the leading contributions to the integral over s' are formed at the region 
M 2 s' ~ M/31/2 ^> 1, provided that M/31/2 <C X- It is interesting to examine an approximation which 
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admits the limiting process E — > 0. For every I < / max in a sufficiently weak electric field, there exists 
such x that M/3? max /2 <C X- Under this condition, we obtain the approximation for finite values of I. 



(O' 



01 



m case 



Note that at low temperatures the terms e Ai<C>/3/ Re jjp and e M<C>/3 ' ReT^J of order e ( M M 
a> 1/2 decrease when / increases. Therefore, any accuracy of calculation for the sums in ([75]) and (|77|) 
can be provided by a finite number of summands. Accordingly, at \^' \ <C M, in order to estimate the 
leading contributions it is sufficient to take into account only the first summand, at I = 1. However, even 
for large values of chemical potentials, \fJ>^ | ~ M, a sufficient accuracy is provided by summation up to 
a corresponding finite number l est . Accordingly, we suppose that l es t < ^ m ax- Thus, the approximated 
expressions for Re^ (t))g and Re(T M „ (t))g at low temperatures have the form 



Re(i M (t))^ = EE(- 1 )' +le * ( ° J ^ 
J* = I ds du h (s, u) Y (s, u) (s, u) ; 







Re(T M , W )^^^(-l) i+1 e^ (0 T^, 
C=± i=i 

rx/M 2 poo 

T* v = I ds du h(s,u)Y (s,u)t^ (s,it) , 
Jo Jo 



(95) 



where h (s, u), Y (s, u), j u (s, u), and t^ L (s, u) are defined in (|8Tj) and ([82)1 . 

Let us find the leading contributions depending on the electric field for expressions (|95|) . These quan- 
tities depend on two dimensionless parameters that characterize the influence of electric field: \qE\ /M 2 is 
a parameter determined by the strength of constant electric field, whereas qE (t + T/2) /M is a parameter 
depending on time t and determined by the increment in the kinetic momentum of a particle in electric 
field by the time instant t (i.e., the increment that accumulates during the time interval t — t\= t + T/2). 
The first parameter \qE\ /M 2 is small in a weak electric field. We shall expand all the functions in the 
integrals (|95|) as series in the powers of this parameter and examine only the leading terms of the expan- 
sion. In addition, the value \qE\ (t + T/2) /M is not necessarily small for a sufficiently large time interval, 
t + T/2. Thus, the dependence on this value in the following expression is taken into account precisely. 
For simplicity, we assume that the magnetic field is weak, so we can neglect the biquadratic terms of the 
order q 2 \EB\ /M 4 . Therefore, with accuracy up to the leading term depending on the electric field, we 
obtain 



TO — - 
[11/ 



oc 

ds 
o Jo 



du h°Y w j™ , 



ds / du 0° + Sh) Y w (*° t + St^) 
o Jo 



(96) 



Here, we have neglected the exponentially small contribution and extended the range of integration over 
s to infinity, +oo; besides, 



-u/4 



2\fnu 



p , Sh 



-u/4 



2 y /nu Air 2 



■ exp 



-MV 



As 



D. 



D = 

T~W = 

<^33 : 
^00 : 



ifiiy 



w_ 

12s 3s 2 

\2 



ifiiy 



- + NT 

s 



jf = -6lq*E(t + T/2)T W /s, Sh 



, Y w = exp 



(qE) 2 (t + T/2) 2 t w u 



u + (j3iy/s 



r3 „2 



u+((3l) 2 /s : " 
(qE) 2 ( s t 2 w s > 



(qE) 



2 ^3 
W 



6s 2 



+ (qEY (t + T/2) 
s \ s 



stkk+stu, st M 



'(^) 

1 M 2 (q E ) 2 ( S -T W f 



fc=l,2,3 



(97) 



where t° and p are given by 
on temperature. Let us do it now. 



Note that in expressions (|9T|) we have not imposed any conditions 
Namely, since at low temperatures the leading contributions to the 



1Q 



integral over s are formed at the region of s 
the powers of u, 



(31/ (2M), the function tw can be expanded as a series in 



TW 



1 



su 



su 



{pi) 



(98) 



We consider only the leading contribution to the expansion in inverse powers of (3M. To calculate this 
contribution, it is sufficient to take into account the first three terms of expansion (|98|k moreover, in 
the coefficients at time-dependent functions it is sufficient to take into account only the first term of the 
expansion. Therefore, at low temperatures the function Y w in the leading approximation has the form 



Y w = Y WL =exp 



WL 



(qE) 2 (t + T/2) 2 



{fit) 



(99) 



Using expansion (j9"5j) and expression (|9l)j) . we can calculate the integrals over u in ([M]) . To this end, 
it is convenient to use the generating function $o (a), defined by the integral 



$ {a) 



1 



2^ Jo 



du 



-w(l/4+a) 



-1/2 



Calculating this integral according to formula (|139|) . we obtain <i>o (ct) = (1 + 4a) 
derivatives with respect to the parameter a, we obtain a formula we need for calculation: 



(100) 



Then, taking 



$n (a) = 



1 



2^ 



du _ e ~u(l/i+ a ) u n = ^_^n <9™$0 (a) 



.1/2 



da r ' 



(101) 



In the case under consideration, we have a = (qE) 2 (t + T/2) 2 (jfij ■ 

Until the T-constant field has switched on, the system in question is at the state of thermal equilibrium. 
When the T-constant field switches on at the instant t\ the particles at the initial state begin to accelerate. 
The corresponding current density and EMT change with time until the moment £2, when the T-constant 
field turns off. Nevertheless, such a non-equilibrium state during a certain interval of time until the 
increment of kinetic momentum is sufficiently small, \qE\ (t + T/2) /M <C 1, is not very different from 
the initial state of thermal equilibrium, and it makes sense to compare these states. It is exactly this 
condition that guarantees the existence of a limiting process E — * 0. Therefore, it is important to 
examine such a case in more detail. On the other hand, it is of special interest to examine the case of 
asymptotically large values of time, when (t + T/2) /M 3> 1. In this case, the presence of electric 
field manifests itself to the maximum. In what follows, we examine two limiting cases. 

Small increment of kinetic momentum When the increment of kinetic momentum is small, \qE\ (t + T/2) /M <C 
1, the kernels in (|96p can be expanded in its powers. To this end, we take into account the fact that 
a = implies that the function (a) is expressed in terms of binomial coefficients, $0 (0) = 1 and 
&n (0) = 2" • 1 • 3 • 5 • . . . • (2n — 1), for n > 1. Using (|96|) after integrating over u, and considering only 
the leading time-dependent contributions, we obtain 



• J u ^ J, 



/>oo 

/ dsp , = T M °„ + AT M „ , 
Jo 



(102) 



where T° and p Q are independent of electric field and given by {SU; besides 



f = -5lq 2 E{t + T/2). 



(103) 



9f1 



In the above representation of T* , the part that depends on the electric field is separated as AT^ U and 
has the form 



4^2 / dsex P 



4s 



At 



Atu 
Ai 33 
At 00 



At 22 = 



Mr 

2s 



Do--\- 



1\ (qE) 2 (t + T/2) 2 



s (piy 



= At u + 



(qE) 2 | (qE) 2 (t + T/2) 2 
2(pi) 2 s 2 



(qE) 2 M 2 D 



(qE) 2 (t + T/2) 2 



(piy 



Do = -■= 



1 (pi) 2 M 2 



MV 



6 12s 



-24 



s 2 

M 2 s 3 



2M 2 



(pi) 2 (pif 



(104) 



Integration in (|104j) leads to the modified Bessel functions of the third kind K v (z), according to formula 
(|9"Tj) . For simplicity, we assume that the chemical potentials are such that e ~( M ^t J - { °)l 3 <§; l. We then 
use the asymptotic expansion 7.13.1.(7) [34] : 



1/2 



1 + \V 



1\ 1 



4/ 2z 



v 2 - 



4 I 8z 



(105) 



As a result of the asymptotic expansion at low temperatures, we obtain the leading contributions of AT^ 
in the form 



AT U = AT 22 = -Q 



M/31 



M (t + T/2) 2 



pi 



AT, 



■s.\ 



AT, 



(10 



3 , A M(t + T/2Y 



2MPI 



31 



235 2 M(t + T/2Y 
256 + /3Z 



1 + 
1 + 



l + O 



Mpi 



( 1 



\MPI 



Mpl 



(qE) 2 t tt \ 



1/2 



-M/3Z 



4tt 2 \2MplJ 

This makes it possible to present the expression for EMT Re(T /J!/ (t))g in (|95[) as follows: 



(106) 



Re(r M „ (*))g = (t, v (t)y e \ E=0 + E E ^ <clA V, 

C=± 1=1 



(107) 



where we have explicitly separated the part (T^ (t))g\ E _ Q , independent of the electric field, which is 
defined in (|55]) . while the terms AT^ depending on the electric field are given by (|106[) . 

With the help of the expression for J* in (|102p . the current density Re(j M (t))g in (p5"]) . with the 
corresponding accuracy, can be expressed in the form 



3 g 2 £(i + T/2) 



M 



J2 n 



(0 



(108) 



where the initial particle density n^' is given by (|88p . Considering only the leading asymptotic contri- 
butions, we obtain as follows: 



i<0 = 2 



M 
2^/3 



3/2 Z 



^ e /3/(^ (C) -A/) 



Z=l 



1 gB/3Z 
W^2M C ° W 



(109) 



If the chemical potential is small, \fJ-^\ <C M, the leading contribution in the sums (|107[> . (|109|) yields 
one term, at Z = 1. In this case, the leading magnetic field-dependent thermal contributions (|M|) in 
Re(T M „ (i))g are rf\ = T 22 = ~ T 33 < 0> t m) > 0> an d au the time-independent contributions for electric 
field-dependent terms in (|107|) are negative. 
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Large increment of kinetic momentum When the increment of kinetic momentum is large, \qE\ (t + T/2) /M 3> 
1, the parameter a in (|100[) is large as well. Then, we calculate the integrals over u in by using the 
asymptotics of $ (a), $ ( a ) — a ~ 1 ^ 2 /2, and obtain 



jx — J f 1 ' 



7~"X = 



2\qE\ (t + T/2) 
(31 

2\qE\ (t + T/2) 



/ dspo s_1 > 
Jo 



where we have used only the first term of expansion (|98[) . At low temperatures, the main contribution 
to the integral over s is formed at the region s ~ (31/ (2M), and thus the leading contributions in (|110[) 
can be found immediately: 

T - - |^| (t + T/2) / **&+A»), (HI) 

where 

a°tn - <* t 22 = , <5°< o = <5°*33 = (qE) 2 (t + T/2) 2 . 
Consequently, using expressions we obtain from representation (|95|) that 



Re<&(*)>3 = -^Msgn^)^ 

C=± 

M 

|g£| (tTT/2) 

Re(T 00 (t)) c e = Re(T ss (t)) c e = \qE\(t + T/2)J2n (0 , (U2) 



Re<Tii(*)>§ = Re(T 22 (t))g = ^ , Re(T n (t))g[ g=0 



where n^' is the initial particle density, given by (| 1 09|) . and Re(Tn |^ = q is the initial value of 
Re(Tu (t))1, defined by flggj) at f < fi. 

4.2.4 High temperatures 

Let us find the leading field-dependent contributions for Re(j M (t))g and Re(T Miy (t))g at high temper- 
atures. In the presence of external magnetic and electric fields, the conditions of the applicability of 
high-temperature asymptotic include the following conditions for the field strengths: 



M[3 < 1 , y/\qB\/3 < 1 , y/\qE\P < 1 . (113) 

In the case of fields that are not very strong, \qB\ /M 2 < 1, \qE\ /M 2 < 1, these subsidiary conditions are 
not relevant; however, they become necessary in the case of strong fields, \qB\ /M 2 > 1, \qE\ /M 2 > 1. 
There are no other necessary restrictions for the field strengths. Under these restrictions, it is convenient 
to use the representation ([75]) . At high temperatures, the leading contributions into the expressions for 

Re and Re T^J are given by the terms J° and T^ v , respectively. These terms emerge from the integrals 
(|80p . In turn, at high temperatures, the leading contributions to these integrals emerge from small values 

of s' , namely, s' ~ ((31) . Therefore, in calculating the leading contributions to all the field-dependent 

functions one can assume that \qB\ s' -C 1 and \qE\ s' < 1. This allows one to use, in the integrals ([50]) . 

expansions in the powers of \qB\ s' and \qE \ s', and then to extend the range of integration over s' to 

+oo, while neglecting the exponentially small contribution. As a result, one arrives at the conclusion 

that and T" , are approximated by the expressions J* and T* in ([TO)) , respectively. Note that at 

high temperatures the leading contribution into the sums (f76|) and (|77p is formed by a finite number of 

summands, when I belongs to the range 1 < I < ( max . A choice for the number Z max is determined by 

-1 



a given accuracy of calculation, and by the fact that Z max <C ( y \qE\(3j and l max -C yy\qB\(3j 

The terms (|96)) are well-decreasing functions of the number I, and therefore in the high-temperature 
limit one can assume that / max — > oo. Thus, the sum is taken over positive integers, as in the case of 
thermal equilibrium, which happens in the absence of electric field. As a result, we obtain the following 
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asymptotic expressions for Re(j M and Re(T MI/ 

oo 

Be(f,(t))? = ££(-iy +1 e» ( °J*, 

OO 

Re(T^(t))§ = ^^(-l) ,+1 e» (0 ^, (114) 

where J* and T*,, are given by expressions (|96[) and (|97p . As usual, in these asymptotics, we suppose 
|m (c) |/3< 1. 

As a whole, we are interested in the leading electric-field-dependent contributions into expressions 
(Hill , as well as in the limiting process -E — > 0, and in the possibility to compare the leading electric- 
ficld-dcpcndent term with the leading magnetic-field-dependent term. Therefore, we assume for simplicity 
that the magnetic field is not stronger than the electric field, \B/E\ < 1; then, within the given accuracy 
of calculation, we can neglect the contributions that contain the magnetic field at powers higher than the 
second, as well as the contributions depending on the product B 2 E 2 . 

At high temperatures, the dimensionless parameter that depends on time (t) and characterizes the 
influence of electric field has the form qE (t + T/2) (3. As in the low-temperature case, this parameter is 
determined by the relative increment of the kinetic momentum of a particle in electric field reached by 
the time instant t; however, the scale on which one estimates this increment is now the temperature 9 
rather than the mass M. 

When the T-constant field switches on, the particles found at the initial state begin to accelerate. 
Nevertheless, such a non-equilibrium state is not very different from the initial state of thermal equilibrium 
when \qE\ (t + T/2) (3-^1, and it makes sense to compare these states. It is exactly this condition that 
guarantees the existence of the limiting process E — > 0. It is clear that at any finite value \qE\ (t + T/2) 
the given condition takes place when [3 — > 0, i.e., it is exactly in this case that we deal with an asymptotic 
expression at high temperatures. Therefore, this case is the most interesting one. Let us note, in 
addition, that this condition implies comparatively weaker restrictions for the value of kinetic momenta, 
\qE\ (t + T/2) <C 1//3 than the values at low temperatures when one requires that \qE\ (t + T/2) <C M. 
For sufficiently high temperatures, ultrarelativistic kinetic momenta can be regarded as sufficiently small. 
The opposite case of asymptotically large values of time at large but finite j3, when \qE\ (t + T/2) (3 ^> 1, is 
a complicated case, when, according to conditions (|1 13[) . the temperature is sufficiently high in comparison 
with M, B, and E; however, at the same time, the temperature is low in comparison with the kinetic 
momentum qE (t + T/2). This case is also interesting, since the presence of electric field manifests itself 
to the maximum. In what follows, we examine these two limiting cases. 



Small increment of kinetic momentum When the increment of kinetic momentum is small, \qE\ (t + T/2) (3 -C 
1, the function Y w in ^ can be expanded as 

Y w = 1 - 5Y , 5Y= (qE) 2 (t + T/2) 2 (r w - s" 1 ^) . 
Let us examine the following approximation: 
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where T° and h° are expressions being independent of electric field; T° is given by ((84)1 : h° and tw are 
given by (|97|) : j° is given by (|103|) : besides 
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where t° is given by (|84)l . St^ and D are given by (|97|) . respectively. 

To calculate the integrals in (| 1 1 5[> and , we need the following functions of a real argument p: 



r , s 1 f°° du u/A 1 

2>V n Jo u / [u+p) 



OR 



By using formula (3.363.2) in [40] . the function I\ (p) can be expressed in terms of the incomplete gamma- 
function: 



(117) 



Differentiating this relation with respect to the variable p, one can find the corresponding representation 
for the remaining functions I n (p) at n > 2. In particular, we have 



(118) 



In the case under consideration, p = ((31) / s. Let us express J* and Atkk for k = 1, 2, 3 in terms of the 
functions /„ (p): 
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where the relation <I>o (0) = 1 for the integral (|100p has been used, and we have neglected the terms 
smaller that the remaining ones when s ~ ((31) . 



Using the representation 9.1(2) in [34] . 

1 p 
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1/2 



du 

2' 4 J ~ 2 J 1 u^J 2 ' 
and formula (j9"Tj) , one can calculate the following integrals: 
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ds, N>2. 



Let us first examine the case n = 1. Calculating the integral over s, we obtain 

2 JV ~ 1 Af 2JV ~ 3 f 00 

./M/3Z 
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z 1 - N K N _ 1 (z)dz. 



(120) 



This representation allows one to extract the leading contribution for M(3 — > 0. The function K^-i (z) 
decreases sufficiently fast at large z, and it is easy to see that the leading contribution is formed at the 
region of small z, z < 5 Z , where 6 Z is a small value, such that M(3l <C 5 Z <C 1. It is sufficient to take into 
account such / that M(3l <C 1, since the contributions \i into the sum (I114| due to the higher values of I 
can be neglected^ Neglecting the small contribution given by the region of z > 5 Z , we now represent Xi 
in the form 



Xl = 



z 1 - N K N _ 1 (z)dz. 



pi J M/31 

This integral can be calculated due to expansion 8.446 [40 for ivjv-i (z) at z <C 1, 

1 /2\ JV " 1 
K N ^(z) = -i-\ (N-2)\. 



(121) 



6 As an alternative, one can sum over I the contribution \i by using the representation H120I I and the summation formula 
with the function K n (xl) |32| . This approach is more effective when one needs to obtain high-temperature expansions 
beyond the leading terms. 
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Then, we finally obtain 
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Xi = " jjt^ ■ (122) 

(27V -3) 03Z) 

in the leading approximation. Expressing I 2 (p) through I\ (p) by using relation (|118p and following the 
same procedure as in the case of calculating X\i we obtain the leading term for \ 2 a ^ N > 3, namely, 

2 27V-3 f(iV-3)! 
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2N-2 
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In a similar way, we can estimate the leading contribution for % 3 . For our purposes, it is sufficient to 

estimate the orders of X3: X3 ~ {ftl)~ A In (Mj3l) at N = 3, and X3 ~ (0l)~ 6 at = 4. Using these results 
for a calculation of the expressions defined by (|116[) and (|119|) . and leaving only the leading contributions, 
we find 
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Proceeding in the same way, we find the leading contribution to AToo as follows: 

AT 00 = ]T AT fefc . 

fc=l,2,3 



(124) 



(125) 



Taking the sum of the terms (| 124[) , in accordance with the representation (|114[) , we use the alternating 
zcta-function 



i(-) = E(-i) w y 



1 



(=1 



and the summation formula (see, e.g., the Appendix of [32]) 

00 
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in which we have fixed only the leading term, in accordance with our approximation. Since r/ (2) = 7r 2 /12, 
we obtain the final result as follows: 
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Even though we examine the case of a small kinetic-momentum increment, the time-dependent terms 
in tL (i) (|127[) may be larger that the logarithmic contributions, because they are limited only by 

the condition that (t + T/2) 2 /[3 2 <C 1/ (qE(3 2 ) 2 . The derivative of these terms with respect to time t 
determines the rate of change for the components Re(T M; , (t))$. 

Let us find the relation of the resulting expression for the current density Re(j M (t))g with respect 
to the initial particle density n'®. The leading terms of the effective Lagrangian A£g\ B=Q (f9"2"]l can be 
expressed as 

2 



it p i — , L 



where formula (|12ip has been used. Then, using (|88p we obtain 
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where 7? (3) = 0, 9015426774, and we have 
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C=± 



Let us compare the time-independent terms in (|127p with the magnetic field-dependent terms in 
{T^iv (t))e\E=Q Using formula (|126p . we find the leading field-dependent term of the effective La- 

grangian A£g (|93]) at high temperatures: 



12tt 2 
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which coincides with previously known results (see, e.g., 
the magnetic field-dependent terms of EMT are given by 



). Then, using formula 



wc find that 
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We can see that the equation of state for the magnetic field in a medium at high temperature is identical 
with the equation of state for this field in vacuum. The increment of the energy density of magnetic field 
caused by vacuum polarization, {Too(t))^. en \ E=0 , and the increment Tq are both negative, which holds 
true in the case of a field of any intensity. On the contrary, the increment of the energy density of electric 
field caused by vacuum polarization, (2oo(t))£ en , is positive in a weak field and is negative in a strong 
one. The time-independent terms of t^„ (t) (|127j) are positive. Therefore, even at the initial period of 
the existence of electric field, when the time-dependent terms of rt v (t) can be neglected, its equation of 
state is essentially different from the vacuum equation and has different forms in the cases of weak and 
strong fields. 

Large increment of kinetic momentum When the increment of kinetic momentum is large, \qE\ (t + T/2) /? 3> 
1, the leading contribution to the integral over u in (II 14[) is formed at the region n> 1. In order to 
find the asymptotic behavior of expressions (jl 14[) . it is sufficient to know a rough approximation of the 
function Y w , namely, 
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Then J* and T* can be approximated as follows: 



JVC 
[IIS 



ds 



ds 



du h°Y j™ 



du (h° + Sh) Y (t° t + <V) 



(129) 



with t° t given by (f54"| . and h°, Sh, St^ given by ([9"T|. where the approximation tw = (fll) 2 /u has been 
used. It is convenient to carry out the integration over u in (|129[) with the help of the generating function 
Yq (a) — exp — (qE) 2 (t + T/2) 2 (pi) 2 a/u , a > 0, which is identical with Yq at a = 1. Using formula 
(|139|) . we obtain 
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This formula, along with the first and second derivative with respect to the parameter a, allows one to 
carry out the integration over u in (|129[) . As a result, we find, with accuracy up to the pre-exponential 
factor, that J* and T* v are of the order exp [— \qE\ (t + T/2) (31]. Consequently, the leading contribution 
is given by the terms J* and T~^ v for 1 = 1. Therefore, considering asymptotically large kinetic momenta, 
one can see that nonzero Re(j M (t))g and Re(T M ^ (t))£ decrease exponentially: 
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A calculation of the pre-exponential factor requires a more accurate approximation and some very tedious 
calculations. 



4.2.5 Strong electric field 



In the limit of a strong electric field, \qE\ ^> M 2 ,(3 2 , \qB\, the leading contributions depending on the 

electric field for the integrals Re j\p and ReT^'J in (|76|) and (|77f emerge from large values of \qE\ s and 
\qE\ t, \qE\ r S> \qE\ s> 1, For simplicity, we assume that B = 0. In addition, we know from Subsec. 
3.2 that the leading time-dependent vacuum contributions, Ke(j fl (t)) p and Ke(T tlu (t)) p do not depend on 
the details of switching on and off when the time interval t — t\ = t + T/2 is sufficiently large, so that 
condition ([38)) holds true. Then, it is instructive to examine thermal contributions for the same large 
interval. In this case, we can use the following approximation: 
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and 
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where 
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and j° is given by expression (|103[) . In this representation, the kernels are products: a function of the 
variable s times a function of the variable m; therefore, these double integrals can be represented as a 
product: an integral over s times an integral over u. Having calculated the integral over s, and using 
formula (f9Tjl for a calculation of the integral over u, we find that 
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where 
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We note that the dimensionless parameter M (t + T/2) is equal to unity for the time period taken 
by light to cross the Compton radius of a particle of mass M. Thus, for a macroscopic period of time, 
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M (t + T/2) 3> 1, one can use the asymptotic expansion (|105p . namely, 
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At low initial temperatures, M/3 3> 1, this asymptotic expansion is valid for any period of time. We 
can see that in the limit of a strong electric field the behavior of the current density Re(j M (t))g and the 
EMT components Re(T M; , (t))g is described by relaxation oscillations. The amplitudes of these oscillations 
are small in comparison with the increasing vacuum contributions Re(j AI (t)) p and Re(T M „(i)) p in (|50p. 
respectively, starting from the instant i when conditions ([38]) are fulfilled. 



5 Discussion and summary 

We have obtained non-perturbative one- loop representations for the mean current density (j M (t)) and 
renormalized EMT (T^ v (t)) ren of a Dirac field in a constant electric-like background. Two cases of initial 
states are considered, the vacuum state and a thermal equilibrium state. In the general case, each of the 
obtained expressions consists of three characteristic terms: 

{?„(*)) = Re(j„(i))S + J*(i), 
(T^(t)) ren = Re(T^(i))£ e „ + Re{T^ (t))g + (t) ; 

see ([oId]) . 

The components Re(T A1J/ (t))J: en given by ([55]) describe the contribution due to vacuum polarization. 
The corresponding components of current density are zero, (j^ (t)) c = 0. These components Re (T MJ/ (t)) c 
can be related with the real part of Heisenberg-Euler Lagrangian. They are local, i.e., they depend on t, 
but do not depend on the history of the process. 

The components Re(j M (t))g and Re(T M1/ (t))g given by f7S]). (1771) . GH|) describe the contribution due 
to the work of the external field on the particles at the initial state. The components J p (t) and t^ v (t) 
given by ([4"8)) . ([50)) . (|69p describe the contribution due to the creation of real particles from vacuum. 
These contributions depend on the time interval t — t\ that passes since the instant the electric field turns 
on, and therefore they are global quantities. 

All these contributions are investigated in detail in different regimes, limits of weak and strong fields 
and low and high temperatures. In all these limiting cases, we have obtained the leading contributions, 
which are given by elementary functions of the basis dimensionless parameters. 

The action of the electric field manifests itself differently for different components of the EMT 
(Tui/(t))ren- F° r instance, in a strong electric field, or in a field which is not strong but acts for a sufficiently 
long time interval, there is the following correspondence: a coincidence of the leading terms for energy den- 
sity, Re(Too (t))g + Tq (i), and the pressure along the direction of the electric field, Re(T 33 (t))$ + t^ 3 (t). 
Note that for vacuum polarization the terms in the mentioned correspondence have the opposite signs: 
Re(Xbo(*))ren = — R e (?33 (*))ren ■ This happens because in the former relation we deal with the equa- 
tion of state for the relativistic fermions, while in the latter relation we deal with the equation of state 
for the electromagnetic field. The dependence of the electric field and its duration for the transversal 
component Re(Tn (t))g + t\ x (t) is quite different from the dependence of the longitudinal component 
Re(?33 (t))e + r 33 (*)• This means that one cannot define a universal functional for the electric field, 
whose variations should produce all the components of the EMT. This is valid not only in case there 
exist particles at the initial state, but also when the initial state is vacuum. Therefore, we can see that 
there cannot be any generalization of the Heisenberg-Euler Lagrangian for the problem of calculating 
the mean values. In our opinion, this fact takes place for any pair-creating background, and is related, in 
such a case, to the existence of nonlocal contributions to the EMT. This explains the failure of numerous 
attempts to examine one-loop effects on the basis of different variants [19l [26] of such a generalization. 

Making a comparison of the constant temperature-dependent components of the mean energy density 
and pressure in a weak electric field immediately after it switches on, we can see that in this case any 
generalization of the Heisenberg-Euler Lagrangian to the case of finite temperature is impossible. In 
other words, even a seemingly small disturbance of thermal equilibrium, produced by a weak electric 
field, exceeds the limits of applicability of the approaches based on thermal equilibrium. 

We have established the restriction \qE\ T 2 -C (which takes place both for the initial vacuum state 
and for a low-temperature initial thermal state) for the strength and duration of electric field under which 
QED with a strong constant electric field remains consistent. Under this restriction, one can neglect the 
back-reaction of particles created by the electric field. On the other hand, there exists another restriction, 



1 -C \qE\ T 2 , see ([55)1 , that allows one to disregard the details of switching the electric field on and off. 
Gathering both restrictions, we obtain the following range of the dimensionless parameter \qE \ T 2 under 
which QED with a strong constant electric field is consistent: 

K \qE\T 2 <7r 2 /2g 2 . 

This inequality is consistent due to it 2 /2q 2 ^> 1, and one can be certain that QED with a T-constant 
field does exist. 

Similar restrictions can be obtained when the initial thermal equilibrium is taken at sufficiently high 
temperatures, f3 \qE\ T <C 1- In this case, we have two inequalities, (3 \qE\ T 3 <C and 1 <C \qE\ T 2 , 
which imply 

1 « \qE\ T 2 « • 
q 2 P\qE\T 

We can see that the upper restriction for \qE\ T 2 is weaker than in the low-temperature case. 

We believe that a similar consideration for an electric-like non-Abelian external field will lead to the 
same restrictions when the created Fermi particles (partons) can be treated as weakly coupled. In the case 
of high temperatures, one can neglect the back-reaction of the created particles on the chromoelectric 
field, in comparison with the contribution from the bosons, which are influenced by temperature in a 
different way (the case of bosons will be examined in another work). 



Appendix 

.1 Separation of particle creation contributions 

Using relations (O, we find 
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(135) 



Taking into account |9|) , we can see that Sg ± describe the contributions due to particle-creation. 
Then, we represent the real-valued terms (j u (t)) e and (T ul , (t)) e in (J57J) as follows: 

(Tfiu (t)) e = Re(T MV (t))§ + Re(T tlI/ (t)) p e , 



(136) 



where 



UAm = iqJ2 teb l >S e ,dx, X >)]\ x=xl , (T^(t)) c e = iJ2 tel^Sg^x')}]^, , 



(137) 



Because of the symmetry j°Sj (x, x'y 7 = 5^ (a/, x), which can be observed in representation (TT9 
we can calculate the real-valued parts of the right-hand sides of (|137j) at X — X 0jS follows: 



tr [luSeA x i x ')]\ x = x , 



lim tr [7 Se,+(x,x')] 
lim tr [j^Se-ix,^)] 



(138) 



and so on. 
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.2 Proper-time representation of S$ 



C 



As has been mentioned, we know the form of the function Sj in (fl"§|) until the moment ti the electric 
field turns on. Therefore, we continue our analysis in the case of t > t\, t' > t\, which we have began in 
the previous section. We are now interested in the components of this functions, which are denoted in 
CGU) as Sex- 

Because of the cutting-off factor, the integrals Sg^ are finite. However, they have derivatives propor- 
tional to the time interval, t—ti, which may be large. Therefore, in the study of conventional contributions 
we have to pay attention to the contributions that grow as t — t\ increases. Note that, as distinct from 
the contributions J v (t) and t p (t), which have been estimated only for sufficiently large intervals t — t\, 
conventional contributions are calculated for any t, < t — t\ < T . 

We remind that the functions Qip n (x) and *ip n (x) are solutions of the Dirac equation in a constant 
field, described in Subsec. 3.1. In the relation ([33)) of that subsection, we have presented the expressions 

for the coefficients G ( - 1 + ) i obtained for the T-constant field, which determine all the coefficients G w f^j 
via the unitarity condition. These expressions for G ( _| + ) at longitudinal momenta p? being large and 



comparable with \qE\ T are different from the asymptotic form at T — > oo. Accordingly, for such large 

values j>3 the coefficients G f cl^) are different from the corresponding asymptotic forms. However, it can 

be shown that under condition (I75|) these differences in expressions for Sg^ can be neglected, and one 

can use, for every p%, the asymptotic form of the coefficients G fcl^)' which holds true for a constant 

field. On the other hand, we need to attract attention to the above-mentioned fact, which will be used 
to obtain the necessary representation, whereas a manifest form of these coefficients in a constant field 
will be unnecessary. 

The energy spectrum of zn-particles at t < t\ has the form 



= e n = \j M 2 +uj+ (tt 3 ) 2 , , tt 3 = +P^j 



where 



n = (pi,n B ,p3,r), uj = \qB\(2n B + 1 - r) , n B =0,1,..., B ^ ; 
n = (p,r), w =p\ +pl, B = Q. 

We expand Nn^ (in) PD|) as a power series of exp j — (3 fe£P — fi^j j and make use of the formula (see 
0D], Eqs. 3.471.9 and 8.469.3) 

r a , i I r du ( u (3 2 l 2 ei\ 
e*p{-flle n } = ^J o ^exp^-^j. (139) 

Since the operator (TL (ti)) 2 = M 2 — (7P^) 2 + (id?, + qET/2) 2 is an integral of motion in the T-constant 
field, the relation 

(H(h)) 2 ±ip n (x) =s 2 n ±'ip n (x) 
holds true at any time instant t. Then, Sg^ (|135[) can be represented as follows: 

S 6 , + {x,x') = l ^T(+)+^ i (a ; )G( + | + ) _1 +^ m (x'), 



Sg^(x,x>) = -i^T(-)_^ n (a:)[G(_r) 



where 



T (C) = £(-l)'e/V°^r-^«p[-H 



du ( u (3 2 l 2 £f 



E\ = M 2 -(jPT) +(-id' 3 + qET/2Y , 
£ 2 = M 2 - ( 7 P±) 2 + (id 3 + qET/2) 2 . 
Representations (1140[) and ([22]) are formally related as follows: 

S e ,± (x, x') = iT^fi^ (x, x') . (141) 



an 



To calculate (j^ (t)) c g and (T^ v (t))g in (|137p according to the prescription (|138[) . we need the expres- 
sions for Sg.+ (x,x') at xq > x' and Se,- (x,x') at xq < x' only. Under this condition, the proper-time 
representation for S^(x,x') has the form 



S T (x,x') = {<yP + M)A*(x,x'), 
±A*(x,x') = A c {x,x'), x -x' ^ 0. 



(142) 



where A c (x,x r ) is defined by the proper-time representations ([55]) and ([56]) : see [TJ [33] [36]. For the 
function Sex { x j x ')i we ml d 



S e , c (z,z') = ( 1 P + M)A c ex (x,x') , 

/>OC 

Ag X (a;,x') = / fex(x,x' ,s)ds , fg X (x,x',s) = r < >' :) f(x,x',s) 



(143) 



where /(x, a;', s) is given by (|56D . 

Thus, we also have proper-time representations for the quantities ( j„ and (T^ The latter 

quantities are finite at finite t — t\. This implies that representation (|143|) can be used in calculating 
the quantities (j^ (t))g and (Tif. (t))g. Then, the quantity (Too (t))g must be calculated via (T^ (t))g as 
follows: 

(Too {t)Y 6 = £ (*)>§ + iM E trft, C (x,i , )| x=ai , , (144) 

fc=l,2,3 C=± 



The latter follows from (|137[) , with allowance for the fact that Sex are solutions of the Dirac equation. 
Then, we represent expressions in (|137j) and (|144[) in the form 

Op (*)>§ = n E tr [7^A^(x, or')] 1^ , 

c=± 



(r 00 («)>§ = E ^ (*)>§ + * m2 E tr Kcfo 1*=* • 

fe=l,2,3 C=± 



(145) 



Following Schwinger, see [8], the kernel f(x,x',s) can be treated as a matrix element of an operator 
in the ^-representation: 



f(x,x',s) = i(x 



exp < — is 



M 



(146) 



Here, \x') is a complete set of eigenvectors of some commuting operators X 1 *, such that X^\x) = x^x), 
{x\x') x — x'). There exist canonically conjugate operators P M and the commutation relations 

—iqF^, where P M „ is the field strength of the external field under consider- 



X",Pu 



i5f,, 



P P 



ation. Using the representation (|146[) . we rewrite fex(x,x', s) (|143|) as follows: 

fexix, x>,s)= 1 £(-l) 1 e^ ( ° -L= ^ ^j- 2 e-^g{x, x> , s, r), 



(147) 



where 



g(x,x',s,T) 



exp 



T (id 3 + qET/2) 2 f(x,x',s,T), T={l3lf/u, 



f(x,x',s,r) = i exp \ — t 



M 2 - 7 P 



exp < — is 



M - 7 P 



(148) 



and 



( 7 p) 2 = ( 7 p||) 2 + ( 7 ^±) 2 , P± = (o,A,A,o) , p,| = (p ,o,o,p 3 ) 



Since the operators ( 7 P|| ) and ( 7 Pj_ ) commute, we present f(x, x', s, r) as follows 



f(x,x',s,T) = i(x 



exp < is I 7 P|| ) — i (s — it) 



M z - 7 P 
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Applying the Schwinger operator techniques [5] to this expression, we find 
f(x, x', s, t) = /|| (a; || , xi , s)f ± (x± ,x' ± ,s~ it) , 



f\\( x h x '\\> s ) = ex P 
fo(x 0) x' ,s) = 



iqA\\ + - (0:3 — x' 3 ) qE coth (qEs) 



fo(x ,x' ,s) 



iqE 



47r sinh (qEs 



exp 



qEs-/ -/ 3 - - (x - x' Q f qE coth (qEs) 



-iqB 



4-7T sin [qB (s — it)] 



f±(x ± ,x'j_, s - it) 
iqA± = -1— (xi - x x ) (x 2 + x 2 ) 



exp 



/ \2 



iqA ± -i(s- it) (M 2 - qBT? 



iqB (xj_ - x' ± ) 
4 tan [qB (s - it)] 

(149) 



The action of the exponential operator in (1148[) is determined by the expression 



ip(s, t) = exp -t (id 3 + qET/2)' 



exp 



Using the integral representation 

r +00 



we find 



% 2 

iqA\\ + - (x 3 — x 3 ) qE coth (qEs) 



<P(s,t) 



a 



a — it 
2 



1/2 



exp < iqA\\ — (x + x + T) 2 



qE\ tci 



2 I a — ii 



b 2 + i2bTqE (x + x' + T) 



(150) 



where 
Therefore, 



4 (a - it) 

a = a (s) = (qE) 1 tanh (qEs) , b — x 3 — x' 3 . 
g(x,x', s,t) = Lp(s,T)f (x Q ,x' , s)f±(x±,x' ± ,s- it) . (151) 
.3 A convenient representation for weak-field approximation 

First of all, we note that at finite temperatures the kernels in (|76|) and (f77|) have singularities due to the 
zeros of (a — it). If we introduce a new variable s', 



then 



s = is + Co , co = (qE) 1 arctan (qEr) 
_ x qE [1 + tan (qEs') qEr] 



(a — ir) 



tanks') l + (qEry 



(152) 



The singular points of (a — ir) 1 are s' = zkirk, k = 0, 1, 2, . . .. Let us deform the integration contour 
over s in integrals (I76|) and (I77|) so that it leaves the origin and proceeds downwards along the imaginary 
axis until the point —ic\, c\ = (it/2 — cq) \qE\ , without touching it, and then proceeds in parallel to 
the axis Res towards the positive infinity. Then, expressions (|76p and (|77l) can be reorganized as follows: 





— J a + J b 




/•ci-0 




= ~* \ 




Jo 








- r°i 




Jo Jo 



ds 



du [h(s,u)Y (s,u)j fl (s,u)]\ s= _ 



du [h(s ) u)Y(s,u)j fl (s,u)]\ 



s—s — i(ci—Q) ' 



(153) 
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and 



rr\a I rpb 

(•ci-0 



T£, = -i ds I du [h{s,u)Y {s,u)t jJlV {s,u)]\ s= _ iS 

Jo Jo 



T° u = / ds du [h(s,u)Y(s,u)t^(s,u)]\ s=s _ i(ci _ 0) , (154) 



^6 _ 

/O "'O 

One can verify that the quantities 



DC: 



poo pO 

5J£ = -i du / ds [h(s,u)Y(s,u)j fl (s,u)]\ s= _ ig , 

Jo Je"(co-0) 
poo pO 

5T^ V = -i du ds[h(s,u)Y(s,u)t IJ ,v(s,u)]\ 3= _ iS , 

Jo Je^(c o -0) 

are imaginary; therefore the real-valued parts of J° and T°„ can be represented as 

Re j; =ReJ';, Re T* v = Re T£ , 

J'° = J« + 8J«, T^=T^ + 5T^. (155) 

We now change the variable s by s' = s + cq in the integrals for J'^ and T'^ v . Then, one can see that 
the contributions from the integration region < u < uq over u in these integrals are imaginary. Then, 
changing the variable u by it', 

u' = u — uq , uo = qE (PI) 2 cot (qEs') , 
and introducing the notation Re — J° , ReT'° v — T® , we represent J® and T" as §EU\i . 
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